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CRITERIA ON CONTRACTIONS FOR ENTROPIC DISCONTINUITIES 
OF SYSTEMS OF CONSERVATION LAWS 

MOON-JIN KANG AND ALEXIS F. VASSEUR 


Abstract. We study the contraction properties (up to shift) for admissible Rankine- 
Hugoniot discontinuities ofnxn systems of conservation laws endowed with a convex en¬ 
tropy. We first generalize the criterion developed in |47], using the spatially inhomogeneous 
pseudo-distance introduced in [50|. Our generalized criterion guarantees the contraction 
property for extremal shocks of a large class of systems, including the Euler system. More¬ 
over, we introduce necessary conditions for contraction, specifically targeted for interme¬ 
diate shocks. As an application, we show that intermediate shocks of the two-dimensional 
isentropic magnetohydrodynamics do not verify any of our contraction properties. We also 
investigate the contraction properties, for contact discontinuities of the Euler system, for 
a certain range of contraction weights. All results do not involve any smallness condition 
on the initial perturbation, nor on the size of the shock. 
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1. Introduction 

In this paper, we develop criteria for the existence of contraction properties of admissible 
Rankine-Hugoniot discontinuities (typically entropic shocks and contact discontinuities) of 
a wide class of systems of n conservation laws endowed with a convex entropy. Consider a 
n x n system of conservation laws 

d t u + d x f(u ) = 0, t > 0, x € M, 
m(0,x) = u 0 (x), 

endowed with a strictly convex entropy rj. 

For the scalar case (n = 1), Kruzkov’s theory [27| shows that the semi-group associated 
with (|l.ll) is contractive for the L 1 norm. For the system case, under small BV perturba¬ 
tion, Bressan, Liu and Yang in (Tl) [36] constructed a L 1 semi-group of solutions. However, 
the L 1 -contraction property does not hold, generally, for systems (See Temple [iS]). 


The Kruzkov’s semi-group is not contractive in L p for p > 1, unless the flux is linear. 
However, Leger showed in m, that any perturbation of an entropic shock wave is con¬ 
tractive in L 2 , up to shift. More precisely, for a strictly convex flux function /, and an 
associated entropic shock ( ui,u r ,a) (i.e., ui > u r ), and for any bounded entropy solution u 
of this scalar conservation law, there exists a Lipschitz shift t H > h(t) such that 



x + h(t)) — S(t, x)\ 2 dx 


is not increasing in time, where S(t,x ) is the traveling wave associated to ( ui,u r ,a) 


( 1 . 2 ) 


S(t, x) 


ui if x < at, 
u r if x > at. 


The shift h{t) depends on the solution u. This contraction property can be extended from 
L 2 , to any relative entropy p(u\S) associated to a convex entropy rj (see next section). Some 
extension to L p , for 1 < p < oo, can be found in Adimurthi, Goshal and Veerappa Gowda 

ra¬ 
in the case of systems (n > 1), this kind of contraction property has been studied in 
jl7l [50] , In [17] the authors developed, in the case of systems, a criterion for contraction of 
admissible (Rankine-Hugoniot) discontinuities. Their criterion is satisfied, for instance, by 
the Keyhtz-Kranzer system with a rotationally symmetric flux. However, it is not applicable 
to many cases, including the Euler system. In [50j, the notion of contraction was extended 
to a family of non-homogenous pseudo-norms, defined for a fixed a > 0, as 


(1.3) 


d(u{t , x), S(t, x)) 


r](u(t,x)\ui) if x < at, 
ar](u(t,x)\u r ) if x > at, 


where S is the traveling wave associated to the studied entropic shock. Notice that the case 
of a = 1 corresponds to the case studied in m • In [50], it is shown, that any extremal shock 
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(i.e. 1 -shock or n-shock) verifies a contraction property, up to a shift, for such pseudo-norms 
with suitable weights a > 0. This pseudo-distance (11,31) (determined by the weight a) does 
not depend on on the solution u. It depends only on the system and the traveling wave S. 

The purpose of this article is to generalize the criterion developed in m to the spatially 
inhomogeneous pseudo-distance introduced in [50] . We first apply our generalized criterion 
to the case of extremal shocks. Then, we develop criteria specific for intermediate shocks, 
and intermediate contact discontinuities. We present two applications of those criteria. 

First, we show that intermediate shocks of the two-dimensional isentropic MHD (which 
is a 4 X 4 system), do not verify the contraction property, for any weight a > 0. For inviscid 
and viscous stability issues for the MHD, we refer to mmmm- 

For the contact discontinuities of the Euler system, it is shown in [46], that the contraction 
property holds for the specific value a = 6 r /0i, the ratio of temperatures on the right, and 
on the left of the contact discontinuity. We show that this cannot hold for a large range of 
other weights a. 

Our criteria depend only on the structure of the system dm and the fixed admissible 
discontinuity. Contrary to the analysis in [[50], it does not involve the study of every solutions 
u, nor the construction of the shift. This simplifies a lot its applicability. The main difficulty 
of the analysis of [50], for the contraction, is due to the construction of the suitable shift 
and weights a. 

The theory of contraction, based on the relative entropy, is valid for large perturbation, 
without smallness conditions. We consider any bounded entropy weak solutions u to mm 
verifying a BVi oc property. The BVi oc property is stronger than the strong trace property 
used in [32., 150]. All our results still hold under the assumption of the strong trace property 
instead of BVi oc . However, following m , we restrict ourselves to the BVi oc case to simplify 
the exposition. Note that the existence theory of entropy weak solutions to the system (ED 
(when n > 3) for large data is open. Strong trace property in the case of scalar conservation 
laws have been widely studied [T5l [ 20 l (28 i IM] 13Ti 1301 [52]. However, in the case of systems, 
the validity of the strong trace property is mainly an open problem. This has been shown 
only for the particular case of isentropic gas dynamics with 7 = 3, for traces in time, in 

E5U 

Our analysis is based on the relative entropy method. It has been first used by Dafermos 
m and DiPerna [21] to show the weak-strong uniqueness and stability of Lipschitz regular 
solution to conservation laws. (See also [TT] [18]) We refer to [16] [25] [31] [32] 03] 06] for 
applications of the relative entropy method to the stability of large perturbation in various 
contexts. This method is also an important tool in the study of asymptotic limits to 
conservation laws. Applications of the relative entropy method in this context began with 
the work of Yau [55] and have been studied in various context (See for instance [ 2 ] [3, 7\ j 8 ] 

m m m m eh m)- 

The rest of the paper is organized as follows. In Section [3] we present the criteria for the 
contraction of admissible discontinuities, and then identify the necessary conditions for the 
contraction. As an application of our criteria, in Section0] we prove the contraction property 
of extremal shocks. In Section [5] using the necessary condition developed in section [3] we 
construct two kinds of criteria preventing the contraction property for intermediate entropic 
shocks, and then present an application to MHD. It turns out that one of the criteria still 
hold true for intermediate contact discontinuities. In Section [ 6 ] as an application of this 
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criterion to gas dynamics, we find a range of weights, for which there is no contraction for 
2-contact discontinuities of the full Euler system. 

2. Preliminaries 

In this section, we present our framework, and basic concepts and properties, needed for 
the analysis in the following sections. 

2.1. General framework. We consider a n x n system of conservation laws: 

(2.4) dtu -I- d x f(u) = 0, t > 0, x € R, 

which is endowed with a strictly convex entropy 77 , thus the system is hyperbolic on the 
state space where 77 exists. Here, the flux /, the entropy 77 and associated entropy flux q 
are assumed to be all defined on an open convex state space V C R ra and of class C 2 (V), 
and the following compatibility relation holds on V: 

n 

djq = diVdjfi, 1 <j<n. 

1=1 

which is conventionally rewritten as 

Vq = VryV/, 

where the matrix V/ denotes (d 


As already mentioned, we have in mind the application of our criteria to the gas dynamics. 
For this reason, we need to extend the phase space V to a suitable subset of the boundary 
of V, to handle the points corresponding to vacuum states. Thus, we introduce as in [53] : 

U = {it € R n | 3uk € V, lim it*. = if, limsup 77 ( 11 *,) < 00 }, 

k—too A:—»0O 

and extend the entropy functional 77 on U by 

77 ( 12 ) = lim inf 77 ( 11 ). 

WEV, U^U 

In 3 x 3 full Euler system, V = (0, 00 ) x R x (0, 00 ) denotes a set of non-vacuum states of 
density, momentum and energy, while U = VU {(0,0,0)} includes the vacuum state (0, 0,0). 
In general case, U is still convex and 77 is convex on U (See [53]). We here restrict our 
study to bounded entropy solutions 11 to (12.41) . whose values are in a convex bounded subset 
Uk C U, on which the functions /, rj and q are continuous. 


2.2. Relative entropy. For the strictly convex entropy 77 of (12.41) . we define the relative 
entropy function by 

rj(u\v) = 77 ( 11 ) — 77 ( 1 ?) — Vrj(v) ■ (11 — v), 

for any u € U ,v € V. 

Since 77 is convex on U and strictly convex in V, we have (see m) 

77 ( 11 ! 7 ;) > 0, u£U,v£V, 


and 

77 (it|i;) = 0 4 => 11 = v. 

Thus, the relative entropy 77 (it|i;) is positive-definite and convex in the first variable it. 
However it looses the symmetry unless 77(it) = |it| 2 . Nevertheless the relative entropy is 
comparable to the square of L 2 distance on any bounded subset of U as follows. 
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Lemma 2.1. For any bounded set B <ZU and compact set fic V, there exists C\. C '2 > 0 
depending on B and 12 such that for any u € B with v € 12, 

C\\u — v \ 2 < rj(u\v) < C^lu — v\ 2 . 

The proof of this lemma can be found in [ 32] 153] . Notice that this lemma also holds for 
all (u,v) G 12 2 , for any compact set 12 e V. 


As mentioned in Introduction, we are interested in studying the contraction property 
of a bounded entropy weak solution as any perturbation of admissible Rankine-Hugoniot 
discontinuity. 

We say that u is an entropy (weak) solution of (11.11) if a weak solution u to (jl.ip satisfies 
the entropy inequality 

(2.5) d t r)(u) + d x q(u ) < 0, 

in the sense of distributions. On the other hand, the equality above holds when u is a 
Lipschitz solution to (11.11) . 


For given ui 7 - u r , we say that ( ui,u r ,a) is an admissible Rankine-Hugoniot discontinuity 
if there exists cr € R such that 


f(u r ) - f{ui) = a(u r - ui) 
q(u r ) - q(ui) < a(rj(u r ) - r]{ui)). 


Equivalently, this means that the discontinuous function u(t,x ) defined by 


u(t, x ) 

is an entropy weak solution to ( 11 . 11 ) . 


ui if x < at, 
u r if x > at, 


For any constant vector v € R n , if u is an entropic weak solution of m, then rj{u\v) is 
a solution in the sense of distributions to 

(2.7) d t ri(u\v) + d x q(u, v) < 0, 

where q{u, v ) is the relative entropy flux defined by 

q{u,v) = q(u) - q(v) - Vrj(v) ■ (/(it) - f(v)). 

This can be derived directly from (11.11) and (12.51) . 


2.3. Spatially inhomogeneous pseudo-distance. For a given weight a > 0, using the 
relative entropy, we consider the pseudo-distance d by 


d(u(t, x), S(t, x)) 


rj(u(t,x)\ui ) if x < at, 
arj(u(t,x)\u r ) if x > at, 


where S(t,x ) denotes the fixed shock ( m,u r ,a ). This pseudo-distance is spatially inhomo¬ 
geneous for a / 1. Based on this pseudo-distance, it has been shown in [50] that there 
exists suitable weight a > 0 such that contraction of extremal shocks holds up to Lipschitz 
shift a(t) in the spatially inhomogeneous pseudo-distance: 

/ CO 

d(u(t, x + a(t)), S(t, x))dx, 

-CO 
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which is equal to 


/ a(t)-\-at /*oo 

rj(u(t,x)\ui)dx + a / r/(u(t,x)\u r )dx. 

-oo J a(t)-\-at 

From now on, we study the contraction properties by using this pseudo-distance denoted 
by (12.81) as 

/ h(t) /*oo 

rj(u(t, x)\ui)dx + a / r](u(t,x)\u r )dx. 

-oo J hit) 


This pseudo-distance (12.81) (determined by the weight a) does not depend on Uk ■ That 
is, it does not depend on any quantitative property of bounded entropy weak solution 
perturbed from admissible discontinuities. On the other hand, the shift h(t) depends on 
the perturbation and is estimated by 


\h'(t)\ < C K , \h(t) - at\ < C K Vt\\u 0 - S\\ L 2, 


where Ck is a constant depending on Uk- 


3. Entropy criteria for contractions 

In this section, we present a general theory for contraction of admissible discontinuity 
for any characteristic fields. First of all, we generalize the criteria developed in m via the 
spatially inhomogeneous pseudo-distance (12.81) . We then give the necessary condition for 
the contraction property in Theorem 13.21 Following the heuristic observation of Serre and 
Vasseur in m, we consider the generalized conditions for contraction as follows. 

Definition 3.1. For a positive constant a > 0, we say that an entropic Rankine-Hugoniot 
discontinuity ( ui,u r ) is relative entropy stable with respect to weight a (in short, a-RES) if 
( ui,u r ) satisfies the following entropy conditions: 

• ((HI) : For any u in T, a := {u \ rj(u\ui ) = ar](u\u r )}, 

D sm {uiy, u ) := aq(u, u r ) - q(u, uf) < 0. 

• (H2) : For any entropic discontinuity (u-,u+) of speed cr± satisfying 

r](u-\ui ) < arj(u-\u r ) and r/(u+\ui ) > ar](u+\u r ), 

Drh{ui, t ,u±) := aq(u + ,u r ) - q(u-,ui) - a±(ar](u + \u r ) - rj(u-\ui)) < 0. 

As a variant of a-RES, we say that an entropic Rankine-Hugoniot discontinuity ( ui,u r ) is 
strongly relative entropy stable with respect to weight a (in short, a-SRES) if (ui,u r ) satisfies 
(HI) and a slightly stronger condition (H2*) than (H2) as follows: 

• (H2*) : For any entropic discontinuity (u_,u+) of speed a± such that 

q(u-\ui ) — ar](u-\u r ) and arj(u+\u r ) — q(u+\ui ) have the same sign , 
Drh(ui t r, u±) < 0. 

Remark 3.1. 1. For the unit weight a = 1, the meaning of 1-SRES is the exactly same 
as one of the terminology ‘RES’ used in m- Thus, the above definition of a-SRES is a 
generalization of ‘RES’. 

2. The condition that rj(u-\ui ) — ar](u-\u r ) and aij(u + \u r ) — g(u + \ui) have the same sign 
has the same meaning as that U- and u+ are separated by the (n — 1)-dimensional surface 
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S a . Note that the convexity and C 2 -regularity of r] implies the convexity and C 2 -regularity 

of 

R := {u | r](u\ui) < ar](u\u r )}, 
thus C 2 -regularity of surface £ a = dR. 

Definition 3.2. For a positive constant a > 0, we say that an entropic Rankine-Hugoniot 
discontinuity ( ui,u r ,a ) satisfies a-contraction if for any bounded convex subset Uk ofU, 
there exists a Constance Ck (depending on Uk) such that the following holds true. For any 
entropy weak solution u € Uk H BVi oc ((0,T ) x R) n (with possibly T = oo) of the system 
m with initial data uq satisfying E(uq, 0) < oo, there exists a Lipschitz function h(t ) with 
h( 0) = 0 such that the pseudo distance E a (u(t), h(t)) is non increasing in time, i.e., 

(3.9) E a (u(t), h(t)) < E a (u(s),h(s)), a.e. t > s > 0, 
moreover, for every 0 < t < T, 

(3.10) \h'(t)\<C K , \h(t) - at\ < C K Vt\\u 0 - S\\ L 2 , 

where S(x) = ui for x < 0 and S(x) = u r for x > 0. 

Remark 3.2. As already mentioned, the weight a > 0 defines the pseudo distance E a (u(t ), h(t)) 
for the contraction, which does not depend on any quantitative property of the perturbation, 
i.e., Uk, but only depends on the system (ED and the fixed discontinuity ( ui,u r ,a). On 
the other hand, the boundedness of the set Uk is for the control of Lipschitz shift h(t ) such 
as (13.101) . Indeed by Lemma UTli E a {u(t),h(t )) is equivalent to \\u(t, ■ + h(t)) — S\\ L 2 , thus 
the contraction property implies 

\\u(t,- + h(t))-S\\ L 2 <C K \\u 0 -S\\ L 2 . 

This induces the estimate (13.101) . Its proof can be found in [50]. Therefore, we will not 
mention about (13.101) in the sequel. 

3.1. The a-SRES implies the a-contraction. We here show that for a given weight a > 0, 
the a-SRES is a sufficient condition for the a-contraction. This is a generalization of the 
main theorem (corresponding to a = 1 ) in [ITT] . 

Theorem 3.1. For a > 0, if the entropic discontinuity ( ui,u r ) is a-SRES, then ( ui,u r ) 
satisfies a-contraction. 

Proof. The proof is almost same as that of the main theorem in S3- For the reader’s 
convenience, we give a variant of the proof of Theorem 2.1 in S3- For given e > 0, we 
define a function V e : R n —> R by 

f [aq(u\u r ) - q(u\ui) - e}+ 

(3-11) V £ {u) = < ar](u\u r ) — r]{u\ui) 

{ 0 if u € S a , 

where X a = {u \ rj(u\ui) = arj(u\u r )}. 

Since ( ui,u r ) is a-SRES, V £ is Lipschitz on R n . Indeed, that is true by the continuity of q 
and the definition of a-SRES. More precisely, since aq(u\u r ) — q(u\ui) — e < — e for u € T, a , 

V £ = 0 on a neighborhood of £ a . 
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We now consider an entropic weak solution u € L°°((0, oo) x R) n fl BVi oc (( 0, oo) x R) n with 
E(uq,0) < oo. Then we define an approximated curve h £ as a solution to the ODE 


(3.12) 


h £ {t) = V £ (u(t,h £ (t))), 

M 0) = 0, 


in the Filippov sense [22]. 

For solvability of (13.121) . we have the following lemma as in [ 37j : 


Lemma 3.1. There exists a Lipschitz solution h £ to (|3.12l) such that 


(3.13) 


H^elll/ 00 < ||Fe |]l°° j 

h £ (t) € I(v £ (u-), V e (u+)), a.e.t> 0, 


where u± := u(t, h £ (t)±), and I(a, b) denotes the interval with endpoints a and b. Moreover, 
if U- / u+, the (u-,u+,h e ) is an admissible entropic discontinuity, that is, 


(3.14) 


/0+) - f( u -) = h £ (u+ - U-) 

q(u + ) - q(u-) < h £ {rj(u + ) - rj(u-)), a.e. t > 0. 


The proof of this lemma is the exactly same as that in B3. because its proof only need 
the regularity of u and V £ , but not the definition of V e itself. Indeed in [31] 32 s [37]) it is 
shown that ()3.12l) has a solution satisfying ([3.131) . We refer to [37] for details of its proof. 
We use (12.71) to derive that the entropic weak solution u € BVi oc (( 0, oo) x R) n satisfies 

d r h c(t) r°° 

—E a (u(t),h £ {t)) = / dtrj(u(t,x)\ui)dx + a d t ri(u(t, x)\u r )dx 

at J —oo Jh e {t) 

+ h £ (t ) (rj(u(t, h(t)-)\ui) - ar)(u(t, /i(t)+)|u r )^ 

/ h e (t) ro o 

d x q(u(t, x)\ui)dx — a / d x q(u(t,x)\u r )dx 
-oo J h £ (t) 

+ h £ (t) (rj(u(t, h(t)-)\ui ) - arj(u(t, /i(t)+)|u r )^ 

< aq(u + \u r ) - q(u-\ui) - h £ {t)(ap(u + \u r ) - rj(u-\ui)) 

=-D £ , 

where u± := u(t, h e (t)±). Let us show that D e < 0. 

For a.e. t > 0 such that U- = u + , by ([3.131) . we have 

h £ (t) = V £ (u±), 

which implies together with (13.111) that 

D e = aq{u+\u r ) - q{u-\u t ) - V £ {u±)(ap{u + \u r ) - rj(u-\ui)) 

= aq{u+\u r ) - q{u-\ui ) - [aq(u + \u r ) - q{u-\ui ) - e} + 

< e. 


On the other hand, for a.e. t > 0 such that U- j - there are two cases as follows: 

i) rj(u-\ui) — arj(u-\u r ) and arj(u+\u r ) — rj(u+\ui ) have the same sign, 

ii) arj(u-\u r ) — rj(u-\ui) and arj(u + \u r ) — iq{u + \ui) have the same sign. 
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Concerning the first case i), we use the fact that ( ui,u r ) satisfies (' H2 *) and (u-,u + ,h £ ) is 
the entropic discontinuity by (13.141) . then it follows from (712*) that 

D e = Drh(ui jT ,u±) < 0, 

where note that a± = h £ . 

For the second case ii), since arj(u-\u r ) — rj(u-\ui ), ar)(u+\u r ) — rj(u+\ui), V e (uJ) and 
V £ (u+) have the same sign, thus h £ € I(V £ (u-),V £ (u + )) also has the same sign. This 
implies together with ()3.14j) that for both v = u- and v = u+, 

D e < aq(v\u r ) - q(v\ui) - h £ (ar](v\u r ) - r](v\ui)) 

= aq(v\u r ) - q(v\u{) - \h £ \\ar](v\u r ) - r}(v\u{)\. 

If we consider v satisfying |I4(n)| = inf(|T4(u_)|, 114 (u+) |), since \h £ \ > |14(u)| by (13.131) . 
we have 

D e < aq(v\u r ) - q(v\u t ) - \V e (v)\\ar](v\ur) - rj(v\ui)\ 

= aq(v\u r ) - q(v\ui) - V £ (v)(aq(v\u r ) - r](v\ui )) 

= aq(v\u r ) - q(v\ui) - [aq(v\u r ) - q(v\u t ) - e]+ 

< £. 

Therefore, it follows from the estimates above that for a.e. t > s > 0, 

E a (u(t), h £ (t)) < E a (u(s), h £ (s )) + (t- s)e. 

Since HV^IIl 00 is uniformly bounded with respect to e thanks to u G L°°((0,oo) x M) n and 
(13.111) . by (13.131) . h £ is also uniformly bounded with respect to s. Thus up to a subsequence, 
h £ uniformly converges to a Lipschitz function h. Hence we conclude that a.e. t > s > 0, 

E a (u(t),h(t )) < E a (u(s),h(s)). 

□ 

3.2. The a-contraction implies the a-RES. The following theorem says that the a-RES 
is the necessary condition for the a-contraction. 

Theorem 3.2. For a > 0, if the entropic discontinuity (ui,u r ) satisfies the a-contraction, 
then ( ui,u r ) is the a-RES discontinuity. 

Proof. Suppose that the entropic discontinuity ( ui,u r ) is not a-RES, by the definition of 
a-RES, at least one of (Til) and (7i2) dose not hold. That is, we assume that one of the 
following conditions holds. 

• (~ Til) : 3 u such that r](u\ui) = aq(u\u r ) and D sm (uiy,u) > 0. 

• (~ TL2) : 3 entropic discontinuity (u-,u + ,a±) such that 

rj(u-\ui) < arj(u-\u r ), r](u + \ui) > ar](u + \u r ) and D RH (ui jr ,u±) > 0. 

We may show that both cases above provide a contradiction with the contractivity (13.91) . 

• Case of (~ Til) : Consider a smooth initial data uq defined by 


[ u 

if x € (-R, R) for some R > 0 

Uo(x) = < Ui 

if x € (—oo, —2R), 

{ U r 

if x G (2 R, oo), 
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where it is a constant vector appeared in (~ Til). 

Then, i? a (rto,0) < oo and the system (11.11) admits the smooth solution u for some small 
time t < Tq, which satisfies 

d rh(t) roo 

—E a (u(t),h(t)) = / d t r](u(t,x)\ui)dx + a d t r](u(t, x)\u r )dx 
at J -oo Jh(t) 

+ h(t) h(t))\ui) - ar](u(t, /i(t))|u r )). 

Since u is smooth for such short time, u satisfies the entropy equality 

d t rj(u) + d x q(u) = 0. 

Thus, for any constant vector v € M n , the smooth solution u verifies 
(3.15) dtr](u\v) + d x q{u\v) = 0. 

We use (13.151) to get 

d f h ® f°° 

—E a (u(t),h(t)) = - d x q(u(t,x)\ui)dx - a / d x q(u(t, x)\u r )dx 
J -oo Jh(t) 

+ h(t)(rj(u(t, h(t))\ui) - arj(u{t , h(t))\u r )^ 

= aq(u(t,h(t))\u r ) - q(u(t,h(t))\ui) 


+ h(t) (v(u(t, h(t))\ui) - ar)(u(t, /i(t))|u r )) . 

For any Lipschitz curve h(t) with h( 0) = 0, since we can choose Tq small enough such that 


u(t, h(t )) = u for all t € [0, To], 
it follows from (~ Til) that for all t S [0, To], 

E a (u{t ), h(t)) = aq(u\u r ) - q{u\ui) + h(t) (ji(u\ui) - aq(u\u r )^ 
= aq(u\u r ) - q(u\ui) 

= D sm (ui t r j u) ^ 0, 


which contradicts with (|3.9p . 

• Case of (~ TL2) : For the entropic discontinuity (u-,u+,cr±) in (~ TL2), we consider 
a initial data no that is discontinuous at x = 0 and smooth on (—oo,0) and (0,oo), and 
satisfies that for some R > 0, 


u 0 (x) 


U- if x € (-R, 0), 
u. |_ if x € (0, R), 
ui if x € (— oo, — 2R), 
u r if x € (2 R, oo), 


In order to freeze the shock speed <r±, we consider a new flux A defined by 


(3.16) 


A(u) = f(u) - cr±u. 


Indeed, by Rankine-Hugoniot condition, we have 

A{u + ) - A(u-) = f(u+) - /(«_) - a±(u + - u-) = 0, 
which means that the speed of entropic shock (u-,u+) to the system 
(3.17) dtu + d x A(u) = 0 
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is zero. Let G be a C 2 -entropy flux such that G' = rf A'. Then, since 

G' = '/(/' - (J ±) = <1 ~ 

we have 


(3.18) 


G = q — a±r/ + C for some constant C. 


Let us consider a weak entropic solution w € L°°((0, oo) x R) n n BVi 0C (( 0, oo) x M) n to the 
system (13.171) with initial data no- Then, there exists small time T* > 0 such that the weak 
entropic solution w is smooth on both (—oo,0) and (0,oo), and has a shock (u_,u+) with 
zero speed at x = 0, and satisfies that for all t < T*, 


(3.19) 


U- if x € (-R/2, 0), 
u .|_ if x € (0, R/ 2), 
ui if x € (—oo, — 3R), 
u r if x € (3 R, oo). 


Let h(t) be any Lipschitz curve with h( 0) = 0. Then, we can choose T* > 0 small enough 
such that 

IM*)l<f t€[0,T.]. 

By the continuity of h, we use the Weierstrass approximation theorem to choose a sequence 
(p n ) of polynomials such that 


(3.20) II Pn - ^||l°°([o,t,]) 0, n-} oo, 

which yields that for sufficiently large N, 

(3.21) \Pn{t)\ < n>N , f€[0,T*]. 

We first show that 


(3.22) E a (w(t),p n (t)) > E a (u 0 ,p n (0)) +tD RH (ui t r,u±), t € (0, T*]. 

Let us begin by noticing the fact that 

R R 

(3.23) w(t,x ) is constant in time t € (0,T*], for x € (— —, —), 

thus w(t,w ) is differentiable with respect to t S (0, T*] for all x € R. This yields 

d rPn(t) roo 

—E a (w(t),p n (t))= d t rj(w(t,x)\ui)dx + a j d t rj(w(t,x)\u r )dx 

J— OO Jp n (t) 

'-v- / >>- v -' 

L I2 

+ Pn(t)(jj(w(t,p n (t)-)\ui) - ar](w(t,p n (t)+)\u r )j . 

- ---V-^ 

h 

(Case of non-constant polynomial p n ) : If the polynomial p n is not constant for some 
n > N, p n has a finite number of zeros as t±, t 2 ■ ■ ■ , tk with 

0 < ti < t2 < ■ ■ ■ < tk < T*. 









12 


KANG AND VASSEUR 


Then p n is either negative or positive on each interval (tj,fj_|_i), where i belongs to one of 
the following classes 

1 < i < k — 1 ifti = 0 and tk = T*, 

0 < i < k — 1 iffi>0 and tk = T* where to := 0, 

1 < i < k ifti=0 and tk < T, * where tk+i '■= T*, 

0 < i < k if 1 1 > 0 and tk < T*. 

First of all, let us consider the case that p n (t) < 0 for all t € (tj,fj + i) for some l. Then the 
solution w is smooth on (0, T*] x (—oo ,p n (t)\. Thus, we apply (13.151) to the system (13.171) 
to get 

rPn(t) 

h = - d x G(w(t, x)\ui)dx. 

J — OO 

Since by (13.191) and (13.21[) . 

(3.24) w(t,p n (t)) = u- for all t € {t h t i+ 1 ), 

we have 

h = -G(u-\ui). 

Since w is discontinuous at x = 0 for 0 < t < T*, we rewrite I 2 as 

. r 


[t r°° 

I 2 = a / dtp(w(t, x)\u r )dx + a / dtr](w(t,x)\u r )dx. 
J Pn (t) 


f r 

4 


By (13.231) . we have 


R 

d t T](w(t,x)\u r ) = 0, (t,x) € (0, T*] x \p n (t), —]. 


Since w is smooth on (0,T*j x [-j,oo) and w(t, -j) = u + for 0 < t < T*, we have 




4 , 4 

P OO 

h = —a / d x G(w(t, x)\u r )dx = aG(u + \u r ). 

J 4 


For I 3 , we use (|3.24f) to get 


h = p n (t){v(u-\ui ) - arj(u-\u r )y 


Thus, we have shown 

E a (w(t),p n (t )) = aG{u + \u r ) - G(u-\ui) + p n (t)(ji(u-\ui) - ari(u-\u r )j, t € 

We here use (|3.18l) to reduce 
aG(u + |tt r ) — G(u_|tt/) 

= a <?(u + ) - o- ± r/(it+) - g(u r ) + o- ± r/(it I .) - drj(u r )(f(u + ) - cr±u + - /(u r ) - cr ± u r ) 
- <?(u_) + a±ri(u-) + q(u{) - a±i)(ui ) + dp{ui){f {uJ) - cr±U- - f(m ) - cr±ui) 

= aq(u + \u r ) - q(u-\ui) - o± ( arj(u + \u r ) - r](u-\ui)j 

= D RH (ui t r,u ±), 
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which provides 


E a (w(t),p n (t )) = D RH (ui t r,u±) + p n (t)(j](u-\ui) - ar](u-\u r )J , t€ {t h ti + 1 ). 

For any r,t G (ti,ti+ 1 ) with r < t, integrating it over [r, t\, we have 

E a (w(t),p n (t )) = E a (w(r ),(t)) + (t - T)D RH (ui >r ,u±) 

+ (Pn(t) - p n {r))(j](u-\ui) - arj(u-\u r )y 

Since E a (w(t),p n (t )) is continuous in time t and p n {ti ) = 0, taking r —>• ti, we have 
E a (w(t),p n {t )) = E a (w(ti),p n (ti )) + (t-ti)D RH (ui tr ,U±) 

+ Pn(t)(ji(u-\ui ) - ar?(u_|xv)). 

Since p n (t ) < 0 for all t G (fz,t;+i), we use the assumption (~ "H2) to get 

(3.25) E a (w(t),p n (t )) > E a (w(ti),p n (ti)) + (t - ti)D RH (ui tr ,u±), te (ti,t t+ 1 ], 

where note that this inequality also holds at t = t] + \ by the time-continuity of E a (w(t),p n (t)). 

On the other hand, for the case where p n (t) > 0 for all t G (t m ,t m + 1 ) for some m, we 
follow the same argument as above. More precisely for all t G (t m ,t m + 1 ), we get 


h = 


R 
' 4 


dtr](w(t,x)\ui)dx + 


rpn(t) 


_R 

4 


d t rj(w(t, x)\ui)dx 


/ 


_R 

4 


— OO 
OO 


<9 x G'(u;(f,x)|itz)Gte = 

I2 =—a I d x G(w(t,x)\u r )dx = aG(u + \u r ), 

J Pn (t) 

h = Pn(t ) ( 17 ( 11 +|uz) - 077 ( 11 + |n r )), 
which provides that for all t G (tmj^m+i)? 

+ Pn(t ) (77(ll+|llz) - 077(11+|ll r )^. 


Since p n (t ) > 0 for all t G (f m ,f m +i), we use the assumption (~ "H2) to get 

(3.26) E a (w(t),p n (t)) ^ E a (w(t m ) ,(t m )) T (t t m ')D R }{(^ui r ,u-i-'), t G (f m ,f m +i]. 

Therefore, combining (I3.25|) and (|3.26l) . we can conclude (|3.22l) . 

(Case of constant polynomial p n ) : If polynomial p n is constant for some n > N, using 
(I3.2ip . we have 


Pn(t) = p n ( 0) G (- j, j) for all t G [0,T*]. 
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First of all, for the case of p n { 0) 7 ^ 0, the claim (13.221) follows directly from the previous 
arguments. If p n ( 0) = 0, we also combine the previous arguments to have 


h = 


f~T f° 

/ d t T](w(t,x)\ui)dx + / d t r](w(t,x)\ui)d: 
J -00 J-% 

d x G(w(t,x)\ui)dx = —G(u-\ui), 


R 
' 4 


' —OO 

. R 


I‘2 = a 


= —a 


f°° 

/ d t i](w(t, x)\u r )dx + a / dtr](w(t,x)\u r )da 

Jo J* 

roo 

/ d x G(w(t, x)\u r )dx = aG(u + \u r ), 


h = 0 , 

which provides (|3.22l) . 

Since tDRH(ui ir ,u±) in (|3.22[) independent of n and (~ V.2) implies 

tDjifj (pi,r, u±) >0 for t > 0 , 
using (13.201) . (|3.22l) and w € L°°([0, 00 ) x M), we show 
(3.27) E a (w(t),h(t )) > E a (u 0 , 0), t € (0,T*]. 

Indeed, since 

|Fl a (u;(t),^(t)) - E a {w(t),p n (t))\ 

fh(t) rpn(t) 

= / p(w(t, x)\ui)dx + 0 / v(w(t, x)\u r )dx 

Jpn(t ) Jh(t) 

< llPn - /l llL°°([0,T*])(ll r ?( u; l^)llL°°([0 ,oo)xR) + «||??(R'|'U r )|| L o o ([o ! oo)xII 
—> 0 as n —> 00 , 


we have 

E a (w(t),p n (t )) -> E a (w(t),h(t)) uniformly in t G [0,T*], 
which implies that for t G (0, T*], 

E a (w(t),h(t )) > -E a (-u 0 ,0) + t -D RH (ui )r: u±) > E a (u 0 ,0). 

Since the Lipschitz shift h is arbitrary, the weak entropic solution u of CLU with uq also 
satisfies 

E a (u(t),h(t)) > E a (u 0 ,0), 0 < t < T*. 

which provides the contradiction with (13.91) . 

□ 


4. On a-CONTRACTIONS FOR EXTREMAL SHOCKS 

In this section, we are going to show that extremal shocks (i.e. 1-shock or n-shock) is 
a-SRES for some a, which implies that they satisfies a-contraction by Theorem 13.11 Even 
though the a-contraction for extremal shocks has been shown by Vasseur in [50], we intend 
to here give an alternate proof as a direct application of the criteria built in Theorem 13.11 
This result goes beyond the known results valid in the class of BV solutions under small 
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perturbation in BV . In the case of small perturbation in L°° n B V, Bressan, Crasta and 
Piccoli in [10] developed a powerful theory of L 1 stability for entropy solution obtained by 
either the Glimm scheme or the wave front-tracking method. The theory also works in some 
cases for small perturbation of large shock (See [331 [9]). On the other hand, Chen, Frid and 
Li in m use the relative entropy to establish the uniqueness and stability of solutions to 
the Riemann problem for the 3x3 Euler system in a large L 1 n L°° n BVi nr perturbation 
(See also mM)- 


4.1. Hypotheses. We suppose the same hypotheses for the system (HD as in [50], which 
are especially applied to the isentropic Euler system and full Euler system. (See [32, [50] ) 

The following hypotheses are related to the 1-shock and other entropic discontinuities. 


• (T~L) : For any fixed ui € V and 1 < i < n, there exists a neighborhood B C V of 
ui such that for any u € B, there is a i-tli Hugoniot curve S l u (s) € U defined on an 
interval [0, s u ) (possibly s u = oo), such that 5^(0) = u and the Rankine-Hugoniot 
condition: 

f(S l u (s)) - f(u) = a % u {s)(S l u {s) - u), 

where a l u {s) is a velocity function. Here, u —>• s u is Lipschitz on U , (s,u) —>• <S*(s) 
and (s,u) —>• cr*(s) are both C 1 on {(s,u) | s € [0, s u ), u € U} : and the following 
conditions are satisfied. 

( 4 - 28 ) (a) < °> ^(0) = A iO)> (b) -^r](u\Sl(s)) > 0, for all s > 0, 

(4.29) X^Slis)) < al(s) < \i(u) < a l u (s), 2 < i < n, s > 0, 
and 

(4.30) ~r^u{s) < 0, (7^(0) = A i(u) for 2 < i < n. 

as 

Regarding the hypotheses for the ra-shock, we just replace (14.28jl and ([4.291) by (14.311) and 

(14.321) as follows. 

• (T~L*) : For any fixed u r € V and 1 < i < n, there exists a neighborhood B C V of 
u r such that for any u € B, there is a z-th Hugoniot curve S^(s) € U defined on an 
interval [0, s u ) (possibly s u = oo), such that S*(0) = u and the Rankine-Hugoniot 
condition: 

f(S l u (s )) - f(u) = a ] 1 (s)(S , ;(s) - u), 

where a l u {s) is a velocity function. Here, u —>• s u is Lipschitz on U , (s,u) —>• ^(s) 
and (s,u) —>• cr^s) are both C 1 on {(s,u) | s € [0, s u ), u € W}, and the following 
conditions are satisfied. 

(4.31) ^ cr™(s) > 0, ^r/(u|S^(s)) > 0, for all s > 0, 


(4.32) A n (S™(s)) > a£(s) > A n (u) > <(s), 1 < i < n - 1, s>0, 

and 

d 


ds 


£7u( s ) < 0, <7^(0) = Aj(it) for 1 < i < n — 1. 
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Remark 4.1. 1. Note that a system verifies the hypotheses ((H) relative to ui € V if 
and only if the system 

d t u - d x f(u) = 0 

verifies (H*) relative to u r € V. In other words, ((H) and (H*) are dual in this way. Thus it 
is enough to show the a-contraction for 1-shocks, because the result of n-shock is obtained by 
applying the case of 1 -shock to u(t,x ) = u(t,—x), which is also entropic solution to ( 11 . 11 ) . 
From now on, we will restrict our arguments to the case of a 1-shock. 

2. The assumptions (a) in (14.281) . (I4.29P and (I4.30P are just due to Liu and Lax entropy 
conditions. The only additional requirement is (b) in (I4.28p . which is a condition on the 
growth of the shock along S^(s) measured through the pseudo-distance. This condition arises 
naturally in the study of admissibility criteria for systems of conservation laws. In partic¬ 
ular, it ensures that Liu admissible shocks are entropic even for moderate to strong shocks. 
Indeed, this fact follows from the important formida : (See also jT7l f30 . 35 ] ) 

Q{Sl{s)) - q(u ) = crl(s)(r)(Sl(s) - ? fiu)) + J ^ ^a l u (t)iq(u\Sl(t))dt. 

In g], Barker, Freistuhler, and Zumbrun showed that the stability (and so the contraction 
as well) fails to hold for the full Euler system if hypothesis (b) in (14.281) is replaced by 

^(^(s)) >0, s > 0. 

This shows that the strength of the shock is better measured by the relative entropy rather 
than the entropy itself. 

3. In [39], Texier and Zumbrun showed that the hypotheses {(H) implies the Lopatinski 
condition of Majda. 

4.2. Structural lemmas. We first present the following structural Lemmas treated in |50j . 
The first lemma provides a kind of triangle inequality for the pseudo metric induced by r?(-1 -) 
and its analogous inequalities, which are useful tools in the following proofs. 


Lemma 4.1. For any u,v,w € W 1 , we have 

(4.33) T]{u\w) + r)(w\v) = r]{u\v) + (V?7(tr) — \/r](v)) ■ (w — u ), 

and 


(4.34) q(u\w ) + q{w\v) = q(u\v) + (V??(u>) — V??(u)) • ( f(w ) — /(«))• 

Therefore, for any a G M, 

q{u\v) — ap{u\v ) = (^q{u\w) — cnri{u\w)^ + ^(rclu) — a r q{w\v) S j 

- (Vr]{w) - V?y(u)) • (f{w) - f{u) - a{w - u)j 


(4.35) 


Proof. The proof follows directly from the definition of the relative entropy and its 

flux < 7 (-|-). Indeed, the following computations hold: 


rj{u\w) + r/{w\v) = (rj{u) — r](w) — '\7r](w) ■ (u — w)^j 
+ (fi]{w) — rj(v) — S7rj(v) ■ (w — v)^ 


= rj{u) — r](v) — Vr}(v) ■ (u — v) + (Vr](w) — S7rj(v)) ■ (w 
= q{u\v) + (' Vr](w ) — Vr/(u)) • (w — u ), 


u) 
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q(u\w) + q(w\v) = ( q(u ) - q(w) - Vq(w) • ( f(u ) - /(»)) 

+ (q(w) - q{v) - Vr?(w) • ( f(w ) - /(u))) 

= ?(«) - g(w) - Vt/(u) • (/(«) - /(«)) + (Vr?(w;) - V^w)) • (/O) - /(«)) 
= g('wl^) + (Vt7 (V) - \/t](v)) • (/O) - /(u)). 

□ 

The following lemma gives an explicit formula concerning the entropy lost at an entropic 
discontinuity (u, S'* (s)) for any i-family. 

Lemma 4.2. For any Rankine-Hugoniot discontinuity (u, S*(s), cr l u (s)) and any vector v, 
we have 

(4.36) q(S z u (s),v) - a t u (s)r](Sl l (s)\v) = q(u, v ) - a l u (s)q(u\v) + J' ^a l u (t)r](u\S l u (t))dt. 
Therefore, for any s > 0, sq > 0, we have 

(4.37) g(5;(s),5*(s 0 ))-<(s)r ? (5) t (s)|5') t (so)) = J* ^l(t)(r](u\Sl(t))-ri(u\Sl(so))^dt. 

In particular, for any u € B as in hypothesis (TL), there exists 6 € (0, ^) and k > 0 such 
that 

.. _ 0 . Q( s l( s ), Sl(so)) ~ °l(s)v(Sl{s)\Sl(s 0 )) < -k\ai(s) - ^(s 0 )| 2 , for \s - s 0 | < 5, 

(4.38) . , . 

q(Sl(s), St(s 0 )) - vl(s)ri(Sl(s)\Sl{s 0 )) < -k\al(s) - <ri(s 0 )|, for |s - s 0 | > 8 . 

The proof of this lemma can be found in [50]. We refer to the work of Lax [30] for the 
estimate (14.361) . And the estimates (14.371) and (14.381) are variations on a crucial lemma of 
DiPerna [21]. Note that the discontinuity (ui, S^(s), a l u (s)) in (14.361) and (14.371) need not be 
extremal family (i.e. 1-family or n-family) from the proof of the relation (14.361) in [50], which 
is obtained directly from the Rankine-Hugoniot condition, moreover, (14.371) is obtained by 
using (14.361) twice. We give the proof of (14.381) in the Appendix for the reader’s convenience. 

For the proof of Theorem 14. 11 we present the following lemma that states both Lemma 5 
and Proposition 2 in [50] as a slightly improved version. In fact, Proposition 2 in [50J says 
that (14.411) holds for all s > 0 satisfying cr u (s) < erg, which means that (14.411) holds for any 
strong shock (u, S^(s), a^(s)) thanks to the assumption (TL) with ([4.281) . It turns out that 
the constraint a u (s) < ao can be removed in the following lemma. 

Lemma 4.3. Let (ui,u r ,a^ r ) be an entropic 1-shock such that u r = S) t; (so), cr^ r = cr^(so) 
for some so > 0, and the corresponding conditions in ([4.281) are satisfied. Then, there exists 
(Jo € (<Ji, r , Xi(ui)), £o > 0, /3 > 0 and a* > 0 verifying the following properties: 

• For any u € B eo (u{), 

cr 0 < Ai(u), 

(4.39) - q(u, ui) + a 0 r](u\ui) < ~/3r](u\ui ), 

q(u, u r ) — Oorj(u\u r ) < — /3r](u\u r ). 

• For any 0 < a < a*, 

(4.40) the ball B eo (ui) contains the convex set R a := {u \ rj(u\u{) < a?7(u|74 r )}, 
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and for any u € R a , 

(4.41) a(q(Su(s),u r ) - crl(s)r](Si(s)\u r )^ - q(u\u t ) + ai(s)rj(u\ui) <0, s > 0. 

Proof. Step A) proof for (|4.39[) : We begin by using (|4.37|) with u = m and s = 0 to get 

q(ui,u r ) - Xi(ui)r](ui\u r ) < 0, 

where we have used the hypothesis (14.281) . Since this inequality is strict, we can choose do 
sufficiently close to Ai (ui) with cr 0 £ {cri,r-> Ai(rq)) and /3 > 0 sufficiently small such that 

q(ui,u r ) - (u 0 - f3)r](ui\u r ) < 0, 

which gives 

q(ui,u r ) - a 0 r](ui\u r ) < -/3r}(ui\u r ). 

Then, using the continuity of q(-,u r ), r)(-\u r ) and Ai(-), we can choose £o sufficiently small 
such that for all u £ B eo (ui), 

q(u,u r ) — aor](u\u r ) < —/3r](u\u r ), and cr 0 < Ai (u). 

On the other hand, we use Taylor expansion at ui to get 

-q(u, ui) + a 0 r](u\ui) = (u - ui) T V 2 r](ui)(aoI - Vf(ui))(u - u t ) + 0 (\u - ui | 3 ). 

Here, since the entropy r] is strictly convex, V 2 r/(rt;) is symmetric and strictly positive, and 
V 2 ? 7 (ui)S 7f(m)) is symmetric. Thus those matrices are diagonalizable in the same basis, 
which gives 

V 2 rj(ui)Vf(ui)) > Ai(nz)V 2 7?(n/). 

This and Ai('iij) > cfq imply that for all u £ B eo (ui), 

—q(u, ui) + a 0 rj(u\ui) < -(Ai (ui) - a o)(u - ui) t V 2 tj(ui)(u - ui) + 0(\u - ui | 3 ) 

< ~(Ai(u/) - a 0 )rj(u\ui) + 0(\u - ui\ 3 ) 

where we have used the smallness of £q in the last inequality. 

Step B) proof for (14.401) : Notice that 

(4.42) 

for a < 1, rj(u\ui) < arj(u\u r ) is equivalent to 

r]{u) < —^—(r](ui) - arj(u r ) - Vrj(ui) ■ u t + aVrj(u r ) ■ u r + (Vr)(ui) - aVrj(u r )) • u), 

1 — a 

where the right hand side of the second inequality above is linear in u. Thus, the convexity 
of i] implies the convexity of R a = {u \ r](u\ui ) < ar](u\u r )}. We take a* < ^ to rewrite the 
second inequality above as 

ri{u\ui) < a ( rj(ui ) - rf(u r ) - Vrj(ui) ■ m + Vrj{u r ) ■ u r + (Vr/(^) - Vr](u r )) • u) 

1 — a 

< Ca( 1 + |n|), for all 0 < a < a*. 

This yields together with Lemma [2~T1 that for all u € R a D B eo (ui), 
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Taking a* sufficiently small as a* < 4A- such that for any a < a* and u G R a n B eo (ui) 


|u - it ;| 2 < ( 7*0 < y, 


which implies that B eo (ui ) strictly contains R a n B eo (ui). Since R a is convex, so connected, 
thus we have 


Ba — Ra n B £q {ui ). 


Therefore, for any a < a*, 

Ba C B £q (ui). 

Step C) proof of (14.411) : First of all, we take oo closer to X\(ui) and Eq smaller than 
those chosen in Step A, such that for all u £ B £Q (u{) 


(4.43) 

O'i(s) < CT 0 , s > y, 

and 


(4.44) 

Ai(o) < 4Ci ki( 2 ) ^i(so)| 


where k > 0 is the constant as in (|4.38[) and C\ is appeared in (14.451) . 

Indeed, (|4.43|) can be justified thanks to the assumption (R), in which ( s,u ) —>• cr u (s) is 
(7 1 -function, and ^cr^(s) < 0 with <r„(0) = Ai(it). 

Let us first show (|4.41|) for all s > We use (J4.35I) to get 


q(Sl(s),u r ) - al{s)rj(Sl(s)\u r ) 

= -(q{u r ,Sl(s 0 )) - al(s)rj(u r \Sl(so))^ + (q(Sfo), S*(s 0 )) - a^(s)r/(S^ (a)|Si(a 0 ))) 
+ (Vrj(ur) - V»7(Si(so))) (/(«r) - /(S'i(a)) - <ri(s)(«r - ^(s))), 


where u r = (sq). By using Rankine-Hugoniot conditions 

f(u r ) ~ f{ui) = a-i,(a 0 )(«r - ui), 
f(Si(s)) - f(u ) = oi(a)(Si(a) - u), 


we rewrite it as 

q(Sl(s), u,.) - crl(s)ri(Sl(s)\u r ) 

= ( - 5 i(a 0 )) + al{s)r](u r \Sl(s 0 ))^ + (q(S*(s), Si(a 0 )) - ^(s)r?( 5 ^(s)|^(s 0 ))) 

+ (VTifa) - Vr?(S’^(s 0 ))) (/(t**) - /(«) - ai(a)(tii - «) + (^(ao) - ffi(s))(«r - uj)) 

= : Ii + / 2 + /3- 


Since u —> a^(s o) and u —>• ^’(sq) are (7 1 and bounded in B £Q (ui ), we have 
k«(ao), 5i(a 0 ))| < C'|S'i I (a 0 ) - ^(s 0 )| 2 < C7|«, - o| 2 , 
? ?(5'i('So)|5'n( s o)) < CIS'^(so) - 5”n( s o)| 2 < C'|«J - o| 2 , 
|V»/(4(ao)) - V77(£i(a 0 ))| < <7|^( So ) - Si(a„)| < C7|«, - «|, 
k^( s o) - criOo)| < C\ui - u\. 
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This yields 

h = -q(Sl t (s 0 ), Si (s 0 )) + (jl(s 0 )v{Sl t (so)I^(sq)) + (<ri(s) - ^(s 0 ))r ? (^ i (s 0 )|S£(a 0 )) 

< C\ui - u\ 2 (l + |<ri(s) - a-i(s 0 )|), 

and 

h = ^(-S^so)) - Vr^S^so))) (/(«0 - /(«) - oi(so)(«i - u) - (<ri(s) - ^(s 0 ))(u; - u) 

+ (ffi,(ao) - oi(so))(«r - «z) - (o-i(a) ~ ^l(s 0 ))(u r - u,)) , 

< C\ui - u\{\ui - u\ + |<ri(a) - cr^(s 0 )||^ - u\ + ^(s) - <ri(ao)|). 

Thus using (I4.d8 p . if |s — so| < <5, we have 

q(Sl(s), Ur) - al(s)rj(Sl(s)\u r ) 

< -k\al(s) - ^(s 0 )| 2 

+ C\u - ui\ 2 {l + |o-i(a) - <7i(so)|) + C\u - uj|ki(s) - ^(s 0 )| 

< C(\u -ui\ 2 + \u- ui \ 4 ) < C\u - ui\ 2 , 

where we have used Young’s inequality and u € B eo (ui). 

If |s - s 0 | > S, then we have 

q(Sl(s),u r ) - o-i(a)^(5i(s)K) 

< -fc|o-i(a) - <ri(s 0 )| 

+ C\u - ui\ 2 (l + |<ji(a) - a-i(so)|) + C\u - «i|ki(s) - ^(s 0 )| 

< C\u - ui\ 2 , 

where we have used smallness of eo as eo < fc. 

Therefore, we use (14.39[) and (I4.43|) to get 

a[q(Sl(s),u r ) - al(s)ri(Sl{s)\u r )) - q(u\u t ) + al(s)r](u\ui ) 

< aC\u - ui | 2 - q(u\ui) + a 0 r](u\ui) 

< aC\u - ui \ 2 + (3r](u\ui). 

Using Lemma 12. II and taking a* small enough such that Ca* < j3 and still R at C B eo (ui ), 
we end up with (14.41(1 for all a < a* and s > 

On the other hand, when s < since 5 < ^ as in (14.381) . we have |s — so| > <L Thus using 
smallness of eo, we have 

q(Sl(s), u r ) - <ri(s)r/(5i(s) \u r ) 

< -fcki(a) - cri(so)| 

+ C\u - ui\ 2 (l + |o-i(s) - o-i(so)l) + C\u - ui\\ai(s) - <ri(s 0 )| 

< -^ki(s) - ffi(so)| + C\u - ui\ 2 . 
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Using (|4.28l) and smallness of £q again, we have that for all s < 4^, 

q{Sl(s),u r ) - crl(s)r)(Sl(s)\u r ) 

< -|ki(y) - o-n( s o)| + C|u - u t \ 2 

< - o-i(»o)|- 

Since u £ R a C B eo (ui), we use (14.28jl and (|4.39l) to estimate 

-q(u\ui) + al(s)r](u\ui) < -q(u\ui) + Xi{u)rj{u\ui) 

= -q(u\ui ) + o- 0 r/(u|uz) + (Ai(tt) - a 0 )q{u\ui) 

(4.45) < (Ai (u) - u 0 )??(tt|^) 

< a(Ai(n) — ao)r](u\u r ) 

< aCi(Xi(u) - a 0 ). 

Therefore, for all s < 4^-, 

0 ( 9 ^(»)> «r) - 0 i(s) 7 ?(‘S«('S)M) - 9(«|«i) + al(s)q{u\ui) 

< ~ “ C ( Al ( U ) “ ^o)) ■ 

Hence we use (J4.44I) to conclude (14.411) . □ 

4.3. Main result. We here revisit the main result in [50j showing that the hypotheses 
(PL) (resp. {PL*)) implies a-SRES , thus a-contraction of 1-shock (resp. n-shock) thanks to 
Theorem 13.11 

Theorem 4.1. Suppose that a system satisfies (fit) and {ui,u r , (Ji, r ) is a 1 -shock such 
that u r = S^ t (so), (Ti r = < 7 * (so) for some so > 0. Then, there exists a small constant 
0 < a* < 1 such that for all 0 < a < a*, (ui,u r ,a^ r ) is a-SRES, thus satisfies a-contraction. 
As a dual result, if we suppose that a system (ED satisfies {PL*) and (ui,u r , a^ r ) is a n-shock 
such that ui = S™ r (so), a^ r = <r” r (so) for some so > 0. Then, there exists a large constant 
a* > 1 such that for all a > a*, (ui,u r ,a^ r ) is a-SRES, thus satisfies a-contraction. 

Proof. By Remark 14.11 we only prove that a given 1-shock ( ui,u r ,oi r ) is a-SRES, which 
implies a-contraction by Theorem 13. II 

• Step A (Verifying {PIT)) : We take <7o, £0 and a* as in Lemma [4.31 such that for any 
0 < a < a*, (14.401) is satisfied, which yields 

E a := {u | r}{u\u{) = arj(u\u r )} C B eo (ui). 

Therefore, using (|4.39l) . we have that for any u € E a , 

D sm {uiy, u) = aq(u\u r ) - q{u\u{) 

= aq{u\u r ) — q{u\u{) + ao{q{u\ui) — aq{u\u r )) 

= a(q{u\u r ) - cr 0 r]{u\u r )j - q{u\u t ) + cr 0 r](u\ui) 

< 0. 
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• Step B (Verifying (712)) : For any 1 < * < n, let (u-,u + ,cr±) be any z-th entropic 
discontinuity satisfying 

(4.46) rj(u-\ui) < arj(u-\u r ) and r/(u + \ui) > arj(u + \u r ). 

Then by Lemma 12.11 the distance |zz_ — ui\ is estimated by weight a as follows: 

(4.47) | U- — ui\ < C\Jrj(u-\ui) < Cy/ar](u-\u r ) < Cyfa. 

For such shocks (zx_, u + , cr± ), we consider 

(4 48) D RH(ui, r ,u±) = a (q(u + \u r ) - a±rj(u + \u r )J + ( - q(u-\ui) + a±rj(u-\ui)J . 

1Z £ 

Case of the 1-shock (u-,u + , a±) : If (zz_, zx + , a±) is a 1-shock, we take £o > 0 and a* > 0 
as in Lemma l4~3l such that (|4.41|) holds for all U- £ R a C B eo (ui), which yields 

DRH(ui,r,u±) < 0. 

where u + = S 3 _(s) and a± = <r„_(s). 

Case of other families i > 2 : If (zz_,zz+, a±) is any z-th entropic discontinuity for z > 2, 
we need to compute C and 7Z in (14.481) further as follows. 

To estimate £, we use Taylor expansion with respect to zz_ at ui, to get 

- q(u-\ui) + Xi(u-)ij(u-\ui) 

(4.49) 

= (u- - ui) T \7 2 q(ui)(\i(ui)I - Vf(ui))(u _ - ui) + 0(\u- - zz z | 3 ). 

Since the entropy rj is strictly convex, V 2 r/(zz;) is symmetric, strictly positive and \7 2 q(ui)\7 f(ui)) 
is symmetric. Thus those matrices are diagonalizable in the same basis, which gives 

X7 2 V (ui)\7f(ui)) > A 1 (zz z )V 2 z ? (u / ). 

Thus together with (|4.47|) . we have 

-q(u^\ui) + \i(u-)rj(u-\u{) < Ca 3/2 , 


which yields 


(4.50) 


C = —q(u-\ui) + \i(u-)rj(u-\ui) + ( a± - Xi(u-))r](u-\ui) 
< (cr± - Xi(u-))rj(u-\ui) + Ca 3/2 . 


On the other hand, we use the identity (14. 35ft and Rankine-Hugoniot condition to have 
TZ = (q(u + \u-) - a±rj(u + \u-)^J + (q(u-\u r ) - (J±rj(u-\u r )^j 
- ( drj(u _) - drj(u r )) • (7(zz_) - f(u+) - a±(u _ - u+)) 

= ( q(u + \u-) — a±r/(u + \u -)) + ( q(u-\u r ) — a±r](u-\u 


= 'h+h- 
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Since (u-,u+,cr±) is entropic discontinuity, it follows from (12.61) that I\ < 0. 

Using (14.351) again, we have 

h = (q(u-\ui) - cr±r)(u-\ui) S ) + (q{ui\u r ) - <J±r){ui\u r )) 

- (drj(ui) - dr)(u r )) • (f{u{) - f(u_) - a±(ui - it_)) 

= : I 21 + I 22 + ^23- 
We use ([4.47)1 to get 

I 21 < C|u_ — ui\ 2 < Ca, 

I 23 < C|u_ — ui\ < Cyfa. 

To estimate I 22 , we apply (I4.37|> with u = ui and s = 0 to the first shock (ui,u r ,cri tr ) so 
that 

q(ui\u r ) - Xi(ui)ri(ui\u r ) = j ' —< 7 ^ (t) [rj{ui (t)) - 77(u/|S , ^(s 0 )))dt. 

By the conditions (f) > 0 and (t)) < 0 in (I4.28p . we have 

(4.51) q(ui,u r ) - Ai(«i)77(u*|u r ) =: -C Jir < 0. 

This and the smoothness of Ai yields 

I 22 = q{ui\u r ) - Xi(ui)rj(ui\u r ) + (Ai (m) - Xi(u-))r}(ui\u r ) + (Ai(a_) - cr±)r](ui\u r ) 

= ~Ci >r + C\u- - ui\ + (Ai(tt_) - cr±)rj{ui\u r ) 

< -Ci >r + Cy/a + (Ai(tt_) - cr±)r](ui\u r ). 

Therefore, by smallness of a, we have 


(4.52) 


77 < -- Ci : r + (Ai(it_) - a±)rj(ui\u r ). 


Therefore we combine (14.501) and (14.521) to get 


(4.53) 


D RH (ui t r,U±) < —~Cl } r + o(Ai(w_) - (T±)r](ui\u r ) 


+ (0"± - Ai(u_))77(u_|^) + Ca 3 / 2 . 

We use (14.47)) to estimate 

7?(a_|u;) < ary(a_|u r ) 

= ar](ui\u r ) + a(? 7 (a_) - rj(ui ) - dr](u r ) • (a_ - a*)) 

= ary(a;|a r ) + aC|a_ — a;| 

< ary(a;|a r ) + Ca 3 / 2 . 

Since Ai(«_) < cr± by (I4.29[) for i > 2, we combine ([4.531) and (14.54)1 to estimate 
Drh(ui !T ,u±) < ~Ci t r + a(Ai(a_) - a±)r}(ui\u r ) 


(4.54) 


+ (a± - Ai(a_))^a7?(a;|a r ) + Ca 3 / 2 ) + Ca 3 / 2 


= -^C z , r + Ca 3 / 2 

< 0 . 
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• Step C (Verifying (712*)) : It remains to consider the case that for 1 < i < n, (u_,u + ,a±) 
is the i-th the entropic discontinuity satisfying 

(4.55) rj(u-\ui) > arj(u-\u r ) and r/(u + \ui) < arj(u + \u r ). 

Then by Lemma 12.11 the distance | u- — ui\ is estimated by weight a as follows. 

(4.56) \u + — ui\ < Cy/r)(u+\u{) < C^/ar](u + \u r ) = Cyfa. 

We apply the identity (14.35(1 to TZ in (14.481) to get 

ft = (q(u+\ui) - <r±v(u+\ u li) + (q(ui\u r ) - cr±r](ui\u r )) 

- (dr](ui) - dr](u r )) • (j(u{) - f(u + ) - a±(ui - u + )j 
= : Ii+I 2 + I3. 

By (I4.56p . we have 

I\ < C\u+ — ui \ 2 < Ca, 

I 3 < C\u+ — ui\ < C\fa. 

We use (I4.5ip . (14.561) the smoothness of Ai to get 

h = q(ui\u r ) - Ai(u;)? 7 (^|w r ) + (Xi(ui) - Ai(u + ))ry(u;|u r ) + (Ai (u+) - o±)r](ui\u r ) 

= -Q t r + C\u + - ui\ + (Ai(tt + ) - a±)rj(ui\ur) 

< -Ci t r + Cyfa + (Ai (u+) - a±)r](ui\u r ). 

Thus by smallness of a, we have 

(4.57) ft < -~Y + (Ai(m+) - cr±)r](ui\u r ). 

On the other hand, applying the identity (|4.35|) to £ in (I4.48p . we have 

£ = ( - q( u _\u + ) + cr±r)(u_\u + )^ + ( - q(u + \ui) + cr±r)(u + \ui)\ 

+ (dv(u+) ~ d v(u-)) ■ (/(«+) - f(u _) - cr±(u + - u-)j 

= : J\ + J2 + J3. 

By (12.61) for entropic discontinuity (u-,u+,a±), we have 

■h <0, J3 = 0, 


thus 

£ < -q(u+\ui) + a±r)(u + \ui). 

We now combine this and (14.571) to have 

DRH(ui >r ,u±) < ~Q >r + a(Ai(u+) - <T±)rj(ui\u r ) - q(u + \ui ) + cr±rj(u + \ui). 

Since (14.551) and (14.561) yield 

V(u+\ui ) < aq(u + \u r ) 

= ar)(ui\u r ) + a(rj(u + ) - rj(ui) - dr)(u r ) ■ (u + - u t )) 

= ar](ui\u r ) + aC\u + — ui\ 

< arj(ui\u r ) + Ca 3//2 , 


(4.58) 
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we use the fact that Ai(it + ) < a± by (I4.29p . to have 

D RH (ui : r,u± ) < ~7;Ci,r + (Ai(u+) - <7±)7?(?z + |u;) + Ca 3/2 - g(u+M + cr ± 7?(it + |^) 
= -|Cz,r - q(u+\ui ) + Ai(«+)r?(u+|uz) + Ca 3 ^ 2 

Moreover, using the same computation as (14.591) . we get 

- q(u + \ui) + \i(u + )t)(u + \ui) 

= (u+ - ui) T 'V 2 r](ui)(\i(u + )I - Vf{ui))(u+ - Ul) + C\u + - ui \ 3 

(4.59) rp r. „ 

< (Ai(n+) - Ai(uz))(u+ - ui ) 1 W 2 rj(u l )(u + - u t ) + C\u + - ui\ A 

< C\u+ - ui\ 3 , 

which yields 

DRH{ui,r, U ±) < ~ 2^’ r + 

< 0 . 


□ 


5. Criteria preventing o-contractions for intermediate entropic shocks 

In this section, as an application of Theorem 13.21 we construct some sufficient conditions, 
which provide no a-contraction for the intermediate entropic shocks. For that, we consider 
the following hypotheses (Mi). 

• (Mi) : We suppose that for some 1 < i < n, the z-th intermediate characteristic 
field (A i,ri) is genuinely nonlinear on U. For a given u € V, suppose that there is a 
C 1 z-tli Hugoniot curve S’ u (s) € U defined on an interval [0, s u ) (possibly s u = oo), 
such that 5^(0) = u and the Rankine-Hugoniot condition: 

f(S l u (s)) - f{u ) = ct u (s)(S , ;(s) - u), 

where cr u (s) is a C 1 velocity function satisfying the Liu entropy condition: 

(5.60) a' u {s) < 0, 
and £T U (0) = A i(u). 

Moreover, we suppose that the system m satisfies the extended Lax condition: 

(5.61) Xi-i(u-) < a < A;_|_i(u_)_), for any Z-th shock (u_, u + , a). 

5.1. Conditions via neighboring genuinely nonlinear fields. The following theorem 
is on the case where there is a genuinely nonlinear field except for the intermediate entropic 
shock. For this case, we construct some sufficient conditions, which provide no a-contraction 
for intermediate entropic shock. 

Theorem 5 . 1 . Consider the system a satisfying (Mi). Let (ui,u r , ai )T ) be a given i-th 
the entropic shock such that u r = S^so), cr U; (so) = cr^ r for some so > 0 and the conditions 
(15.601) and (15.611) are satisfied. Then the following statements holds. 
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• (1) For 0 < a < 1, we assume that there is a C 1 j-th rarefaction curve Ru t (s ) 
with j < i such that Ru t { 0) = ui and the backward curve Ru^(s) of R 3 Ul (s), i.e., 
A j(Ru~(s)) < Xj(ui), intersects with the (n— 1)-dimensional surface S a . Then, the 
entropic shock ( ui,u r ,ai >r ) does not satisfy a-contraction. 

• (2) For a > 1, we assume that there is a C 1 k-th rarefaction curve Ru r (s) with k > i 
such that Ru r { 0) = u r and the forward curve Ru+ (s) of R^ r (s), i.e., A k(Ru+(s)) > 
A k(u r ), intersects with the (n— 1)-dimensional surface S a . Then, the entropic shock 
(■ ui,u ri ai tr ) does not satisfy a-contraction. 

• (3) For a = 1, we assume that one of the assumptions of (1) and (2) is satisfied. 
Then, the entropic shock ( ui,u r ,ai }V ) does not satisfy a-contraction. 


Remark 5.1. The reason of splitting a > 0 into the three ranges as above is motivated from 
the geometric observation. More precisely, the (n — 1)-dimensional surface Si for a = 1 
becomes hyperplane as 

(drj(u r ) - drj(ui )) ■ u = dr)(u r ) • u r - rj(u r ) - dr/(ui) ■ ui + 

which separates ui and u r in the phase space, whereas the surface £ a for 0 < a < 1 (resp. 
a > 1 J is strictly convex that belongs to the region including ui (resp. u r ). Thus, for the 
case of 0 < a < 1 (resp. a > 1), the rarefaction curve Ru t (s) (resp. R u ’+(s)) issued from 
ui (resp. u r ) is likely to intersect with S a . On the other hand, since S a shrinks to ui (resp. 
u r ) as a ^ 0 (resp. a —>• oo), the rarefaction curve Ru(~(s) (resp. Ru+(s)) hardly intersects 
with S a for o» 1 (resp. a <C 1). In fact, the proof of Theorem \6.1\ does not depend on the 
strength of weight a. 

Proof of Theorem 16.11 (1) Case of 0 < a < 1 : Let u be the first intersection point of 
the j-th the backward rarefaction curve Ru^(s) and the surface S a . Let us put u = Rf^(s) 
for some s > 0. We show that the shock (ui,u r ,ai r ) does not satisfies the first condition 
(HI) of a-RES as 

D sm (uiy, u) = aq(u, u r ) - q(u, uf) > 0 . 

For this end, we consider a differentiable function / : M + —>• R defined by 
^ F (s) = a(q(R(l~(s),u r ) - A i (Ri"(s))r ? (R^“(s)|u ? .)) 

- q( F u(f(s),ui) + \j(Rf;~(s))ri(Rl;-(s)\ui). 

Note that we have 

F (s) — D srn (ui )r , u), 

because ar}(u\u r ) = r}(u\ui). 

Since the z-tli shock (ui, S l u (s), er* (s)) satisfies (|5.60p and (j5.61j) . by the inequality (|4.36D 
and the fact A j(u{) < A i(ui), we have 

r s o d 

q(ui,u r ) - (7i, r r](ui\ur) = - J —a l Ul (t)rj(ui\S Ul (t))dt>0, 
m,r > Xi-i(ui) > A j(ui), 


(5.63) 
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which implies the positivity of /(0) by 

F( 0 ) = a(ci(ui,u r ) - Xj(ui)rj(ui\u r )) 

( 5 - 64 ) = a(cj(ui,Ur) - m,r'n{ u l\ u r)) + a(<Jl t r - \j(ui))r}(ui\u r ) 

> 0. 

Let us show 

(5.65) F'(s) >0, 0 < s < s. 

First of all, since -^Ru^(s) = rj(Ru~ (s)), we have 


which yields 


Vf(R{;-(s))±R^(s) = 

F\s) = a[(d V (R^-(s)) - dr]{u r )) (v / (R^~(s)) - A ^“(a))/) ^“(a) 
- ( d V (Ri;-(s )) - dr7(^))(v/«-( S )) - A J ( J Ri’-(s))/)^i-( S ) 
+ ^j(Ru^(s)){ri(Ril~(s)\ui) - ar](R^-(s)\u r )^ 

^i(^ , r( a ))(^(- H «r( a )i“») - a^-wiur)). 


ds 

Here, since 4^\j(R,u i ~( s)) is continuous and A j{Ril^(s)) < Xj{ui) besides 
d_ 
ds 

we have 


-A,«"(s)) = VA.-(^"(a)) • r,-(i2j-(a)) + 0, 


d_ 

ds 


(«)) <0, s > 0. 


Since Ru t (s) = u is the first intersection point of the continuous curve R^ l (s) and the 
surface £ a , we have 

“ av(Ril~(s)\u r ) <0, 0 < s < s. 

Thus we have shown (15.651) . which implies together with (15.641) that 

F’(s) — D sm (ui >r , u) > 0. 

Hence the shock (ui,u r ,ai <r ) is not a-RES, which provides the conclusion by Theorem 13.21 
(2) Case of a > 1 : We follow the same argument as above by considering Ru + instead 
of Ru t . Indeed, we use the same computations to have 
/(0) = -q(u r , ui) + \ k {u r )ri(u r \ui) 

= - (q(ur,ui) - ai, r r](u r \ui) S ) + ( \ k (u r ) - ai, r )ri{u r \ui) 

> 0 , 


/'(») = ^ X j(RuM s ))(jl( R uH S )\ U l) - a V(Ru? OOK)) 


>0, 0 < s < s, 


and 
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where u := Ru+(s) is the first intersection point of the forward curve i?«l + (s) and the 
surface T, a . Thus we have the same conclusion as 

/(s) — D sm {u^ r , u) > 0 . 

(3) Case of a = 1 : Since the proof of two cases above does not depend on the strength 
of o, if one of the assumptions of ( 1 ) and ( 2 ) is satisfied, we end up with no a-contraction 
of the entropic shock (ui, u r , ai )T ). 


5.2. Application to magnetohydrodynamics. As an application of Theorem 15.11 we 
here show that there is no contraction property of certain intermediate shocks for two- 
dimensional (planar) isentropic MHD in Lagrangian coordinates: 


(5.66) 


dtv — d x u = 0 
dt{vB) — /3d x w = 0 
d t u + d x {p + \B 2 ) = 0 
d t w - fid x B = 0, 


where v denotes specific volume, and two-dimensional fluid velocity ( u , w) and magnetic 
field (/3,B) only depend on a single direction e± measured by x. This behavior of two- 
dimensional vector fields with spatially one-dimensional dependence is achieved when the 
initial condition is so. Thus the divergence-free condition of magnetic held of full MHD 
reduces that /3 is constant (See for example [ 6 ] for study on (|5.66p l . As a perfect fluid, the 
pressure p is assumed to satisfies 


(5.67) 


p(v) = v 7 , 7 > 1 . 


This system has an entropy 7 as 

roo 1 2 

V{U) = J p(s)ds + -{u 2 + w 2 ) + 
in terms of the conservative variables U := ( v,q,u,w ) where q := vB. 

For simplicity of computation, we use non-conservative variable W := (v, B,u,w) and 
rewrite (15.661) as a quasilinear form: 

d t W + Ad x W = 0, 


where the 4x4 matrix A is given by 


f 0 

0 

-1 

0 \ 

0 

0 

B 

V 

V 

-c 2 

B 

0 

0 

V 0 


0 

0 / 


where c := \/—p'(v) denotes the sound speed. 

Since the eigenvalues of A solves the characteristic polynomial 

' B 2 + (5 2 , 2 \ x2 , /3 2 


A 4 - 


+ & 


A 2 + 


-c z = 0 , 


we have four eigenvalues 


Ai — —v / o+, A2 — — \Jat-, A3 — y / a_, A4 — ,/a+, 




where 


1 \B 2 + (3 2 


2 


B 2 + (3 2 


2 
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By a straightforward computation, we have the corresponding eigenvectors 


n 


r 2 


r 3 


r 4 


1 ’- (a+ B C2) ,^ /5 +'- 


B 

(a_ - c 2 ) — 

— I- -—-, V®-! 


- 1 


B 

(«+ ~ c 2 ) 
B 


V&+, 


/3(a+ - c 2 ) ^ 

/ ’ 

/j(«- - C 2 ) ^ 
£ v /aT / ’ 
- c 2 ) \T 
BJaZ ) ’ 
/3(a+ -c 2 ) \T 
BJa^. ) 


Here, we restrict our study to the case of 5 / 0. Using the strict convexity of pressure 
p(v), d\i ■ ri > 0 for all 1 < i < 4, which means that all characteristic fields are genuinely 
non-linear. Indeed since 


a ± -c =- 


B 2 + fd 2 


-c 2 ± 


£ 2 + /3 2 


— c“ 


+ 4- 


B 2 c 2 1 


yields 

(5.68) a- < c 2 < a + , 

thus we have 


dX\ ■ r\ = 




- d v a + + 8bol+ 


— c 


, 2 -, 


£ 


and 


-<7a+ = - 


B 2 + /3 2 


2 1 u 2 


+ P + 


(B 2 +0 2 _2 V B 2 +/3 2 „//x , 2B 2 c 2 , 2B 2 p" 

y—v - c li v 2 - p ) + — i 


V / (a±£- c2 ) 2 + 4 


B 2 c 2 


> 


2 B 2 c 2 I 2B 2 p" 

I n} 


( B 2 +P 2 „ a 


> 0 , 


-c 2 +4 




a+ — c 

77 


B 2 c 2 


.2 r j ( «!±i£_ c 2)i + 2£ 


-h 


(«+ - c 2 ) 


> 


V /( a£ J £ -7 2 + 4^ 

= (a + — c 2 ) > 0. 


2 c! 

V 


v /(b!±£_ c 2) 2 + 4 


B 2 c 2 
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Similarly we have dXi ■ r^ > 0 for i = 2, 3,4. 

Let (Ui, U r , a 2 ) be the 2-shock wave satisfying the Rankine-Hugoniot condition: 

~[u] = C2H, 

-(3[w\ = cr 2 [q\, 


(5.69) 


b] + 




Ivi 


02 [u], 

<?2 H, 


where [/] := f r - fi. 

By Lax condition, dX 2 ■ r 2 > 0 implies that —r 2 {Ui) is a tangent vector at Ui of the 2-shock 
curve Sy issuing from [//. Thus since dv ■ (—r 2 ) < 0 and du ■ (—r 2 ) < 0, we have 


(5.70) 

Since 


[v] < 0 and [u] < 0. 


dB ■ (—r 2 ) 



< 0 if B > 0 
> 0 if B < 0, 


we have that [B] < 0 for Bi > 0, and [B] > 0 for Bi < 0. In particular, we here consider 
the case where the 2-shock wave satisfies 


(5.71) either Bi > B r > 0 or B[ < B r < 0. 

Similarly for a given 3-shock wave (Ui,U r , <73), we have [u] > 0, [u] < 0, and [B] > 0 for 
Bi > 0, and [5] < 0 for Bi < 0. Also, we consider the case where the 3-shock wave satisfies 

(5.72) either B r > Bi > 0 or B r < Bi < 0. 

We are now ready to show that for any a > 0, there is no a-contraction of such interme¬ 
diate shocks as follows. 


Theorem 5.2. Let ( Ui , U r ,a 2 ) be a given 2-shock wave of the system (15.66(1 - (15.671) satisfying 
()5.7ip . Then there is no weight a > 0 such that ( ui,u r ) satisfies a-contraction. Likewise, 
this result holds for a given 3-shock wave (Ui,U r ,afi) satisfying (15.721) . 

Proof. First of all, we show that for any 0 < a < 1, the backward 1-rarefaction wave 
issuing from Ui intersects with the three dimensional surface £ a , i.e., 


E a := {17 | rj{U\Ui) = av(U\U r )}. 


Since dv ■ r\ = 1 > 0, v is strictly monotone along the integral curve of ri, which means 
that the 1-rarefaction wave can be parameterized by v. Moreover since d\± ■ r\ > 0, — ri is 
the tangent vector of the backward 1-rarefaction wave Rjf ~, which implies that v decreases 
along Ry'~. That is, v + < vi for all parameters v + of R]}~ ■ Notice that Ry~ is well-defined 
for all v + € (0, v{\, because —r\(W) is smooth for all W £ (0,00) x (R — {0}) x R 2 . Indeed, 
since 


(5.73) 


a + -c =- 


B 2 + fi 2 


-c 2 + 


B 2 + fi 2 


— c- 


2 R2 c 2 1 

2) -f 4 7— 1 


> 0 , 
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we have 
(5.74) 


dB • (-n) = 


(a+ - c 2 ) 


B 


>0 if B > 0 
<0 if B < 0, 


which implies that B + / 0 along Ry due to Bi / 0, thus —r\(W) is smooth for all 
W € (0, oo) x (R — {0}) x R 2 . 

We now use the fact that 

for a < 1, r](U\Ui) < aiq{U\U r ) is equivalent to 

v(U) < - arjiUr) - Vr?(C/ z ) • U t + aVrj{U r ) ■ U r + (Vry^) - aV V (U r )) • U), 

1 — a 

which is rewritten as 

f°° 1 < 7 2 

(5.75) / p(s)ds + -(u 2 + w 2 ) + j— < ci + c 2 (u + q + u + w), 

Jv 2 2v 

for some constants c\ , C 2 . This implies that 

r](U\Ui) < ar](U\U r ) v > c* and |g| + |u| + |u>| < c* for some constants c*, c* > 0, 

since / 0 °° p(s)ds = +oo, and the positive terms on u , w and q are quadratic in the left-hand 
side of (I5.75p . Therefore there exists 0 < »* < c* such that 


fl( R u t ( v *)\Ui)> ap(R^J i (u*)| U r ), 

which implies that Ry~ intersects with £ a for a < 1 , because Ry is a continuous curve 
issuing from U[ € {U \ r](U\Ui) < ar]{U\U r )}. 

On the other hand, we show that the forward 4-rarefaction wave Ry issuing from U r 
intersects with the surface E a for any a > 1. Since dX 4 ■ r 4 > 0 and dv ■ < 0, 

(5.76) r 4 is the tangent vector of the forward 4-rarefaction wave Ry +, 

and the parameter v + decreases along Ry+ ■ Moreover Ry^ is well-defined for all v + € (0, v r \ 
by the same reason as above. 

Let us consider a continuous function 




-. 1 ,- 


We claim that 
(5.77) 

Using (15.761) . we have 
dF 1 (Rj}+(v+)) 
dv+ 


F a (U) := niUpt) - ar,(U\U r ). 

Fi{Ry+(v*)) < 0 for some u* € ( 0 , v r \. 

dR^(v + ) 


= (S7r)(U r ) - V V m) 


= b] + 


. 2 v 2 - 


+ 


dv+ 
v(a -1 — c 2 ) 
Q+ 


+ N\/a+ - M 


v + j3(a + — c 2 ) 


Q+V&+ 


And (15.691) yields 

dF i( R u+( v +)) , , ,—, r t , v+(a+ - c 2 ) ( f3 2 \ 

- — -= ( 2 h_ v ^)H + ---(1-^=) 


h 


O2y[oi+) M 


12 
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Since (|5.73p yields 


(5.78) 


V^T > V~p'( v +) 



—> oo 


ns u_j_ —y 0-t~, 


it follows from (I5.7UI) that I\ —y —oo as v + —y 0+. 
To control I 2 , we use the condition (15.711) . Since 


dB ■ r 4 


(g + - c 2 ) 
B 


> 0 
< 0 


if B > 0 
if B < 0, 


if Bi> B r > 0, £ 


[B] < 0 and q + > 0. Since a + —>• +00 by (15.781) . we have 
h < 0 for v + <C 1. 


This is also true in 
Thus we have 


the case of Bi < B r < 0 because of f 


> 0 and q + < 0. 


dF(R^(v + )) 

dv + 


—>• —00 


as u_)_ —y 0-t - , 


which implies (15.771) . 

Therefore we conclude that Ry + intersects with £ a for any a > 1, because F a (U r ) > 0 and 
F a (F 4 ^{v,)) < Fi (i?4 + (u»)) < 0 for all a > 1. 


Hence for all a > 0, the 2-shock wave (Ui, U r , 02 ) does not satisfies a-contraction property 
thank to Theorem 15 m 

Similarly we use the same arguments as above to show non-contraction for 3-shock wave 
( Ui,U r ,a 3) satisfying (15.721) . More precisely, we can show that the backward 1-rarefaction 
wave Rl~~ intersects with 

Ui 

S a := {U | r)(U\Ui) = ar)(U\U r )} for any 0 < a < 1, 

and the forward 4-rarefaction wave R^ + intersects with S a for any a > 1. We omit the 
details. □ 

5.3. Conditions via neighboring linearly degenerate fields. The following criterion 
is on the case where there are linearly degenerate fields as neighboring families of the 
intermediate entropic shock. 

Theorem 5.3. Consider the system satisfying (Mi). Let (ui,u r ,ai tr ) be a given i-th 
the entropic shock such that u r = 5* z (so), cr Ul ( sq) = ai r for some so > 0 and the conditions 
(I5.6(jp and (15.611) are satisfied. Then the following statements holds. 


• (1) For 0 < a < 1, we assume that for some j < i, the j-th characteristic field is 
linearly degenerate, thus there exists the j-th Hugoniot curve S J Ul (s) as the integral 
curve of the vector field rj with S 3 Ul { 0) = ui such that the contact discontinuity 
(ui,Si l (s),a J Ul (s)) satisfies the Rankine-Hugoniot condition and 

(5.79) = °’ = 

Moreover if there is si > 0 such that 

v(Sh(si)\ui) > arjiSi^s^lur), 


(5.80) 
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then the entropic shock (ui,u r ,a^ r ) does not satisfy a-contraction. 

• (2) For a > 1, we assume that for some k > i, the k-th characteristic field is 
linearly degenerate, thus there exists the k-th Hugoniot curve Sff r (s) as the integral 
curve of the vector field r^ with S^ r (0) = u r such that the contact discontinuity 
(Su r i s )i u n (7 u r ( s )) satisfies the Rankine-Hugoniot condition and 

(5.81) fal(s) = 0, <4.(0) = \ k (u r ). 

Moreover if there is s 2 > 0 such that 

(5.82) r](S^ r (s2)\ui) < ar](S^ r (s2)\ur), 

then the entropic shock {ui,u r ,a^ r ) does not satisfy a-contraction. 

• (3) For a = 1, we assume that one of the assumptions of (1) and (2) is satisfied. 
Then, the entropic shock {ui,u r ,a^ r ) does not satisfy a-contraction. 


Remark 5.2. For example, the condition (15.801) (resp. (15.821) ) is geometrically satisfied by a 
case that C 1 curve (resp. (s)J transversally intersects with the (n—1)- dimensional 

surface £ a . 

Proof of Theorem 15.31 (1) Case of 0 < a < 1 : Let us put := ui, u + := Su_(si) 
and a± := aL_(s 1), then the contact discontinuity (u_, u + , a±) satisfies 

rj(u-\ui ) = 0 < ar)(u-\u r ) and r){u + \ui) > arj(u + \u r ). 

We are going to show that the shock (ui,u r , a^ r ) does not satisfies the second condition 
(7*2) of a-RES as 

Drh(ui,t', u±) := aq[u + ,u r ) - q(u-,u t ) - a±(ar](u + \u r ) - rj(u-,ui)) > 0. 


In fact, since q(u-,ui ) = rj(u-\ui) = 0 by U- = ui, we have 


D RH (uiy,u±) = a(q(u + ,u r ) - (T±r)(u + \u r )j. 


We use the identity ()4.35|) and the Rankine-Hugoniot condition for the discontinuity (u_, u+, a±), 
to get 

D RH (uiy,u ±) = a q(u + ,u t ) + q{u h u r ) + (dr](ui) - dr](u r )) ■ (f(u + ) - f(ui )) 




+ rj(ui\u r ) + (drj(ui) - drj(u r )) • (u + - u t 


= a 


q(u + ,ui) +q(ui,u r ) - a± 


+ v(ui\u r 


where the fact u_ = ui is used above. 

We decompose the above relation into three parts by 


D RH (uiyu±) = a(q(ui,u r ) - Xj(ui)r](ui\u r )^ + a(Xj(ui ) - a±)ij(ui\u r ) 

+ a(q(u + ,ui ) - a±ri{u + \ui )) 

= : ii + h + h- 














34 


KANG AND VASSEUR 


Following the same argument as (|5.63D by using (|4.36p and the fact A j(ui) < we have 


h = a (q(ui,u r ) - ai >r r](ui\u r )) + (<Ji >r - \j(ui))r](ui\u r ) 

Since A j(ui) = <7^(0) = <r^(si) = o± by (|5.79lh we have 

h = a{Xj(ui) - (T±)r)(ui\u r ) = 0 . 

We use (14.36(1 together with (15.7911 to get 


> 0. 


I 3 = a(q{u + ,u-) - a±r](u + \u-)) =a 


f 

Jo 


d_ 

ds 


v 3 u_(s)v(u-\S J u _(s))ds = 0. 


Therefore we have 

Dmi(ui,r',u ±) > 0 . 

Hence the shock (ui,u r ,cr^ r ) is not a-RES, which provides the conclusion by Theorem 13.21 
(2) Case of a > 1 : We follow the same argument as above by considering instead of 
Sl r By (15.821) . we have 

rj(u-\ui) < ar}(u-\u r ) and rj(u+\ui) > 0 = ar](u + \u r ). 

for the contact discontinuity (u-,u+,a±) where u + := u r , U- := (^2) and <J± : = 

<^ + (s 2 ). 

Since q(u + ,u r ) = r)(u+,u r ) = 0 by u + = u r , we have 

DRH(uiy,u ±) = - q(u-,u t ) + a±rj(u-,u l ). 

Using (14.351) and the Rankine-Hugoniot condition of the discontinuity (u_, u + , <r±), we have 


D RH (uiyu± ) = -q(u-,u+) - q{u r ,ui ) + cr ± ^r/(u_|u+) + r/(u r |u;)^, 

where the fact u + = u r is used above. 

We decompose the above relation into three parts by 

Drh{ui,t\u ±) = ~(q{u r ,ui) - X k (u r )r](u r \ui )) + (< r ± - A fc (u r ))ry(tt r .|^) 

- (?(«-,«+) - cr±77(u_|u + )^ 

= : Ji + J2 + J3. 


By the same argument as the first case, we have 

Ji = - (q(u r , ui) - ai, r r]{u r ,ui)^ + (A fc (tt r ) - a^ T )r]{u r ,ui) > 0. 

Since A k{u r ) = a^ r ( 0) = cr^(s 2 ) = a± by (15.811) . we have J2 = 0. 

We use (14.36(1 together with (1 5.81 (1 to get 

Js= ~L ^ a ^+^ r, ^ u+ ^+^ ds= °- 

Therefore we have 

DRH(uiy,u±) > 0 . 

(3) Case of a = 1 : Since the proof of two cases above does not depend on the strength 
of a, if one of the assumptions of (1) and (2) is satisfied, we end up with no a-contraction 
of the entropic shock (ui,u r ,ai >r ). 
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6. Criteria preventing o-contractions for intermediate contact 

DISCONTINUITIES 

In this section, as an application of Theorem 13.21 we consider some sufficient conditions, 
which provides no a-contraction for intermediate contact discontinuities. Then we apply 
the sufficient condition to the full Euler system in order to find a range of weights a, on 
which that intermediate contact discontinuities do not satisfy a-contraction. 

6.1. Condition preventing a- contract ions for intermediate contact discontinu¬ 
ities. In this part, we consider the case that the system fll.ll) has a genuinely nonlinear 
field, which has same conditions as Theorem 15.11 In fact, the proof of Theorem 15.11 still 
holds in the case of the contact discontinuity (ui,u r , a^ r ) instead of shock. Thus we state 
the following theorem without proof. 


• (.4.2) : Suppose that for some 1 < i < n, the z-th characteristic field (Aj, r^) is linearly 
degenerate. Then, for a given z-th contact discontinuity {ui,u r ,a^ r ), there exists 
the z-th Hugoniot curve S l Ul ( s ) as the integral curve of the vector field zy such that 
5^ (0) = ui and S^so) = u r for some so > 0, and the Rankine-Hugoniot condition: 

- f(ui ) = <(s)(S^(s) - u{). 

where a l Ul (s) is a velocity function satisfying the Liu entropy condition: 

(6.83) ^- a l Ul (s) = 0, 0<s<s 0 , 

and cr^(0) = A i{u{), a z Ul (s 0 ) = a^ r . 

We also suppose the extended Lax condition: 

(6.84) A/_i(?z_) < a < A; + i(zz_)_), for any Z-th shock (zx_, u + , a). 

Theorem 6.1. Consider the system satisfying ( A 2 ) and ()6.84[) . Let (ui,u r , ai, r ) be 

a given i-th intermediate contact discontinuity (i.e., 1 < i < n) such that u r = S l Ul (so), 
a\ H (a'o) = <j i t r for some so > 0. Then the following statements holds. 

• (1) For 0 < a < 1, we assume that there is a C 1 j-th rarefaction curve Ru t (s) 
with j < i such that Ru t ( 0) = ui and the backward curve I?4 _ (s) of R J Ul (s), i.e., 
A j(Ru~(s)) < Xj{u{), intersects with the (n— l)-dimensional surface £ a . Then, the 
contact discontinuity (ui,u r ,a^ r ) does not satisfy a-contraction. 

• (2) For a > 1, we assume that there is a C 1 k-th rarefaction curve Ru r (s ) with k > i 
such that Ru r ( 0) = u r and the forward curve Ru+ (s) of Rff r (s), i.e., \k(Ru+ (s)) > 
Xk(u r ), intersects with the (n — 1)-dimensional surface £ a . Then, the contact dis¬ 
continuity (ui,u r , ai !r ) does not satisfy a-contraction. 

• (3) For a = 1. we assume that one of the assumptions of (1) and (2) is satisfied. 
Then, the contact discontinuity {ui,u r ,(T^ r ) does not satisfy a-contraction. 
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Application to gas dynamics As an application of Theorem 16.11 we here find weights 
a > 0, on which the contact discontinuity of full Euler system for a perfect gas does not 
satisfy a-contraction. The 3x3 full Euler system reads 

( d t p + d x (pv) = 0 

(6.85) < d t (pv) + d x (pv 2 + p) = 0 

[ d t {p{\v 2 + e)) + d x (p(\v 2 + e)v + pv) = 0, 

Here, the equation of state for a perfect gas is given by 

(6.86) p = (7 - 1 )pe, 
where 7 > 1. 

We consider the conservative variables p, q := pv and £ := p(\v 2 + e) and put u := (p, q,£), 
and the entropy 

(6.87) rj(u) = (7 — l)plnp — pine, 
where e is given by e = ^ — in conservative variables. 

By a straightforward computation (See for example [45]), the three characteristic fields are 
given as 

Ai = v - c, n = ( - p, c, -p _1 p) , 

A2 = V, r 2 = ( - d e p, 0, -dpp') , 

A 3 = v + c, r 3 = (p,c,p -1 p) , 
where c := 7/7(7 — l)e denotes the sound speed. 

Thus since d\ ■ r* > 0 for i = 1,3 and d \ 2 ■ r 2 = 0, the first and third characteristic fields 
are genuinely non-linear, whereas the second characteristic field is linearly degenerate. Let 
(ip, u r ) be 2-contact discontinuity, then we have 

(6.88) vi = i v, Pi=p r - 
We refer to [45] for these relations. 


Theorem 6.2. Let ( ui,u r ) be the 2-contact discontinuity for the system (16.851) - (16.861) . 
Then, there is no weight a for the range of either 

e r 

0 < a < —, a > 1, when e/ > e r , or 
(6.89) &l er 

a < 1, a > —, when e; < e r , 

ei 

such that ( ui,u r ) satisfies a-contraction. 

Proof. We show that the assumptions of Theorem 16.II are satisfied for a belonging to (16.891) . 
For that, we separate the range of weight a > 0 into two cases. 

i) Case of either 0 < a < — when e/ > e r or a < 1 when ei < e r . For such a, we 
are going to show that the first backward rarefaction wave R u \ issuing from ui intersects 
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with the surface £ a . As in [45], it turned out that the first backward rarefaction wave Ru\ 
starting from ui can be parametrized by the pressure as follows. 


R},\ is parametrized as pi < p+, 


u t 


(6.90) 


p .|_ = 4 >(p+; pupi ), where <p is increasing in the first argument, 


v+ = vi 


L 


p+ dp 

pi P c ' 


where c is the speed of sound, in particular c = 1 / 7(7 ~ l)e for a perfect gas. That is, along 
Ru\~, the pressure p + and density p + increase while the velocity v + decrease. Then, using 

(6.91) 

we have 


pc = p\J 7(7 — l)e = sj^pp by (16.861) and (16.901) . 

1 r p 

\fl JPi 


V+ = Vi 


r p + dp_ 
V^Jpi Vpp' 


Multiplying by y / p+, we have 

— — 1 f p+ [pT dp 

VP+v+ = VP+vi -— / J - - 

\/7 J pi \ P VP 

, _ 1 [ 

<^ VI j 


\/7 jpi 

rP+ dp 

V^Jpi Vp 


= \/ P+Vl - —( y/p^ - s/pi), 


where we have used p + > p thanks to (|6.90l) . 
If p+ is bounded in p + , then we have 


which implies 
(6.92) 


sfp- py + —> —00 as p + —> 00 , 

either p + —>• 00 or 8 + —>• 00 as p + —>■ 00 . 


This means that the either density or energy should be unbounded along the first backward 
rarefaction wave R.// issuing from ui. 

On the other hand, it turns out that ui is enclosed by the surface £ a in the case of either 
0 < a < I 7 when ei > e r or a < 1 when ei < e r as follows. 

Recall (|4.42l) . i.e., 

for a < 1, p(u\ui) < ap{u\u r ) is equivalent to 

v(u) < —C rj(ui ) — ap(u r ) — Vp(ui) ■ ui + aVp(u r ) ■ u r + (Vp(ui) — aVp(u r )) ■ u ). 

1 — a 


Since the entropy (I6.87P can be rewritten as 


V(u) = 7Plnp - pln(£ - d-), 


(6.93) 
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using pe = £ — we have 


(6.94) 


d p r](u ) = 7 (ln p + 1 ) - ln(pe) - —, 


2 e 


dgV(u) = 

e 

<9 £ p(m) = —. 

e 

We use (16.861) . (| 6 . 88 j) and (|6.94l) to compute 

(Vp(rq) - aVr](u r )) ■ u = ( d p rj(ui ) - ad p p{u r ))p 


(6.95) 


By (|4.42|> . (16.931) and ()6.95f) . we have 


+ (. d q rj(ui ) - ad q rj(u r ))q + 


a 1 


ei 


where 


np - ln(5 - |^) - ci) p + c 3 £ < c 2 q + c 4 , 


1 


ci = 


-( d p r]{ui ) - ad p r](u r )), 


l — a 

c 2 = --( d q rj(ui ) - ad q r](u r )), 


1 — a 

1 / a 1 


C3 = 


1 (z v ey / 

c 4 = 1 ( 7 ?(m) - ap(u r ) - Vrj(ui) ■ ui + aVp(u r ) ■ u r ). 

1 — a 

Since C 3 > 0 in the case of either 0 < a < f 21 when e; > e r or a < 1 when e; < e r , we use 

®l 

2 

C 25 < _)_ C5 p t)y Young’s inequality, to get 


np - In (5 - —) - ci - c 5 ) p + —£ < c 4 . 


c 3 , 


Notice that c 4 is constant while 

L —>■ 00 as £ —>■ 00 or p —>■ 00 , 

which implies that {it | ry(ri|ri;) < ar](u\u r )} is bounded. Thus ui is enclosed by E a . There¬ 
fore by this and (16.921) . we conclude that the first backward rarefaction wave Ru\ has to 
intersect with the (n — l)-dimensional surface £ a . 


ii) Case of either a > 1 when e/ > e r or a > ^ when ei < e r . For such a, we 

ez 

3 

use the same argument as previous to show that the third forward rarefaction wave R v ’ r 
issuing from u r intersects with the surface E a . Since the third forward rarefaction wave 























39 


3 + 

Ru’ r starting from u r can be parametrized by the pressure p + as: 

R ^’ r + is parametrized as p r < p + , 

p + = 4 >(p+; p r ,p r ), where is increasing in the hrst argument, 


(6.96) 


v + = v r + 


I p+ dp 
Jpr P C ’ 


the pressure p + , density p + and the velocity v+ all increase along . Then by (16.911) , 

1 f p+ dp 


v + — v r = 

Multiplying by y / plj7, we have 

VP+ V + = VP+ V r + 


\/7 J Pr VPP' 

1 f p+ I p + dp 


Pr 


> \/~P+Vr + 


V7. 

1 f p+ dp 

\fl Jp r \[P 


p Vp 


= \fP+V r + — (VP+ - VP?)> 

where we have used p + > p thanks to (16.961) . If p + is bounded in p + , then we have 

^y~p^v + —>• oo as p + —>• oo, 

which implies 


(6.97) either p + —> oo or £ + —>• oo as —>• oo. 

This means that the either density or energy should be unbounded along the third forward 
rarefaction wave Ru issuing from u r . 

On the other hand, by the same computation as previous case, u r is enclosed by £ a . Indeed, 
for a > 1, rj(u\ui) > ar}(u\u r ) is equivalent to (|4.42[) . besides C 3 > 0 in the case of either 
a > 1 when e; > e r or a > f 21 when e/ < e r . Therefore by this and (16.971) . we conclude that 

3 1 

the first backward rarefaction wave R Ur has to intersect with £ a . □ 


Remark 6 . 1 . Recently in [36], the authors has showed that the 2-contact discontinuity 
of full Euler system (16.851) satisfies a-contraction for the specific weight a = — in the 
borderline of weights as in (16.891) . More precisely, Euler system is written in Lagrangian 
mass coordinates is considered and it turns out that the 2 -contact discontinuity ( ui,u r ) 
satisfies jj--contraction, where 61 and 9 r denote the temperature corresponding to the left 
end state ui and right one u r , respectively. In fact, since e = c 6 for some positive constant 
c for the perfect gas, the weight a = ^ is equal to a = ff, moreover, since it is known 
that the a-contraction property is invariant under a transformation between the Eulerian 
coordinates and Lagrangian coordinates in [47] (See also [44] and [53] ). Serre’s result also 
holds full Euler system (16.851) in Eulerian coordinates. 


7. Appendix 

We here give the proof of (14.381) in Lemma 14721 Since f^a^s) and ^?/(tt|5„(s)) are both 
continuous and non zero at s = sq by (14.281) . there exists 5 G (0, sq) such that for all s with 
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|s - S 0 | < S, 
(7.98) 


T, ~ l 'iW 


< 


ds 


O'i(so) 


-^v(u\Sl(s)) - ^(u|S£(s 0 )) 


1 

< - 

“ 2 




which can be estimates as 
d 

(7.99) 


<7^(s) =-—(T^(s) 

dt uy ’ dt y 


< 




■^rj{u\ s l(s)) = ^-r](u\Sl(s)) 


> 


1 


ri(u\Si(s 0 )) 


Thus, using (14.371) and (17.991) . we have that for all s with |s — so| < <5, 

q(Sl(s), Sl(s 0 )) - al(s)r](Sl(s)\Sl(so)) = j ^lit) (r]{u\Sl{t)) - r](u\Sl(s 0 )fj dt 


< - 


< -- 


dt 
d 

Jt auK 


^ a u( s o) ^l(u\Si(s 0 ))\s- s 0 \ 2 . 


i(»°) -^v(u\ s U s o)) ( t-s 0 )dt 


Since it follows from (|7.98|) that 


d i/\ 3 d i. . 

d, aM ~ 2 * <r “ (So> 


we have 


ki(s) - <ri(»o)| < § ^ a n( s o) \s - sol, 


which yields that for all s with |s — so| < <5, 

Q(Sl(s), Si(s 0 )) - al(s)r](Sl(s)\Sl(so)) < -fci|<ri(s) - ^(s 0 )| 2 , 
where k\ is a positive constant as 


fci = min 

U^B 


On the other hand, let us show (|4.38l) for |s — so| > S. For all s with s < so — 5, we use 
(14.37P to get 

q(St(s),Sl(s 0 )) - ai(s)r](Si{s)\Si(s 0 )) = ^l(t)(j](u\Sl(s 0 )) - r}(u\Sl(t))) dt 

+ j s - ri(u\Sl(t))jdt 

= : h + h- 


If we consider a positive constant c± satisfying 

(7.100) ci < nun (ji(u\Sl(s 0 )) - r](u\Sl{s 0 - 5))^ , 
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using (14.281) . we have 

h < J ^°l(.t)(ri(u\Sl(so)) - r/(u\Sl(s 0 - S))jdt 

< -ci|ct(s 0 -5) - <r(s)| 

< -ci|ct(s 0 ) - cr(s)| + ci|o-(s 0 ) - o-(s 0 - <5)|. 

Since I 2 < 0 for all u € U, we choose c\ sufficiently small such that (17.1001) and 

ci|o-(s 0 ) - 0 -O 0 - <5)| < - min I 2 , 

u&B 


which yields 

q(Sl(s), Sl( 3o )) ~ °l(s)r]{Sl(s)\Sl(s 0 )) < -ci|a(s 0 ) - cr(s)\. 
Similarly, choosing c 2 > 0 sufficiently small such that 

c 2 < min (rj(u\Sl(s 0 + (5)) - rj(u\Sl( 

rso+5 d 


u£B 


C 2 |o-(s 0 ) - cr(s 0 + <5)1 < - min 

u&B J SQ dt 


(v(u\Sl(t)) - ri(u\Sl(s 0 ))) dt, 


we have that for all s with s > so + 6, 

q(Sl(s), 5i(s 0 )) - al(s)r](Sl(s)\Sl(so)) < -c 2 |a(s 0 ) - <r(s)|. 
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